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Tray Drying Fine 


Drying rates finely powdered aluminum soaps 
and diatomaceous silicas were determined tray 
dryer. The weight solid the tray was auto- 
matically recorded. expected for fine powders, 
plots drying rate against moisture content were 
concave downward. The data were correlated 
drying rate and the free 
moisture content lb. water per lb. dry solid. 
straight line was obtained plotting log for 
aluminum stearate against the reciprocal absolute 
temperature and similar relationship was assumed 
hold for the other materials which were dried 
two temperatures only. Intercepts the temperature 
plots seemed depend upon bulk density average 


particle size but general correlation was possible. 


term drying usually refers the removal liquid 

from solid applying heat and carrying away the 
vapor gas either natural forced convection. 
tray drying, heated air passed over the wet solid 
which placed shallow trays supported fixed 
removable racks. large proportion the air, about 
80-90%, recirculated reduce the amount heat 
required. Tray drying batch process far the 
solid concerned but the air flow continuous. 

the investigation reported this paper, mixtures 
water and fine powders were dried tray suspended 
from balance the top the dryer. The eight 
the tray was automatically recorded and the drying rate 
correlated with moisture content. 


THEORETICAL 


The theoretical and fundamental concepts solids 
drying have been summarized Marshall 
refer more particularly the drying granular solids. 

general, liquid moves diffusion through homo- 
solids such sodium soaps and capillary 
action granular fibrous materials and 
some fine powders such pigments and clay, the mech- 
anism depends upon the moisture content. 

Water held the spaces between the particles 
granular solid three different depend- 
ing upon the amount water present. the 
state, all the pore spaces are completely filled with water 
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the funicular state, the particles are covered 
continuous film water but there some air the 
interstitial spaces. the pendular state, the small 
amount water present held the points contact 
the particles. Evaporation from granular solids during 
the capillary state causes the liquid surface recede 
between the particles the surface and the resulting 
curvature produces suction. evaporation proceeds, 
the suction increases. When exceeds certain value 
(the entry suction air enters the solid through the 
larger capillaries replace the water carried the 
surface through the smaller capillaries and the water 
changes the funicular state. Evaporation the sur- 
face and transfer water capillary action continue 
until the reserves liquid within the solid dry up. When 
the water the surface the pendular state, evapora- 
tion occurs zone that gradually recedes from the 
surface. Water moves the evaporation zone and the 
vapor must diffuse through the solid the surface and 
then the air stream. 

When experimentally determined drying rates are 
plotted against moisture contents, lines resembling curve 
Figure result. Section represents the con- 
stant rate period, the first falling rate period and 
the second falling Point defined 
the critical moisture content and point the equili- 
brium moisture content. 

During the constant rate period, evaporation from 
surface film water that replenished diffusion 
capillary action rapidly depleted. some 
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HEAT TRANSFER 
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Figure drying rate and heat transfer 
cient with moisture content. 
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point will longer possible for water reach the 
surface rapidly can evaporated and the dry- 
ing rate will decrease. This critical moisture content 
depends upon the solid, the particle size, and the bed 
thickness well the properties the air stream. 
The first falling rate period called the period un- 
saturated surface because drying occurs the 
surface but diminishing fraction it. The rate 
drying frequently linear with moisture content 
during this period and depends upon the air stream pro- 
perties, the critical moisture content, and the rate during 
the constant rate period. The period unsaturated 
surface drying may not occur all may occupy 
the entire falling rate period. The second falling rate 
period begins hen the pendular state exists the sur- 
face. The rate drying now depends the move- 
ment water within the solid and independent the 
velocity and humidity the air but, course, still 
affected the air temperature. Although often 
calculations the drying time, had generally been 
accepted that the portion the drying rate curve 
was concave upward for both granular 
solids. Newitt and 13) obtained concave 
upward curves for some granular solids but with fine 
particles microns) the curves resembled N-2 Figure 
Oliver and Newitt concluded that beds fine 
particles, the suction pressures were that the 
continuous threads liquid were broken the release 
air and vaporization liquid. Liquid water, there- 
fore, could not carried the surface capillary 
action but, instead, was evaporated below the surface. 

Several authors have that the overall heat 
transfer coefficient between the air and the solid 
constant long evaporation occurs the surface 
even though the drying rate decreasing but that 
decreases the plane vaporization recedes into the 
solid (curve Figure 1). 


Quantitative relationships 


Heat transferred the solid convection from 
the air, conduction from the sides and bottom the 
tray, and radiation from the tray above, 

the constant rate period, the total heat transferred 
related the drying rate 


The coefficient for the transfer heat flat plate such 


2/3 
cpVp 


for Reynold’s number above 15,000. this case, the Rey- 
nolds number expressed 


Re = (3) 


However, higher coefficients are obtained tray drying, 
probably the rougher surfaces and, for the 
air-water system, Treybal recommended 


for between 500 and 6000 (sq. ft.). The radia- 
tion heat transfer coefficient may estimated from 


The temperature the solid surface not ordinarily 
known but during the constant rate period the 
transfer equation 


may combined with 


arrive the equation 


hr 


This may solved simultaneously with the saturated 
humidity curve give and Hg. The drying rate can 
then calculated. These equations are applied later 
evaporation from free water surface. 

Many authors have presented equations for the falling 
assumed linear relationship between drying rate and 
distribution assumed the start the falling rate 
period, equation (9) approximate form the dif- 
fusion equation. 


this equation, the diffusivity water the solid 


and changes with moisture The em- 
pirical linear equation expressed 

(10) 
where 


According Friedman and Marshall “), equations (10) 
and (11) apply coarse granular solids and certain 
materials high moisture contents while equation (9) 
applies single phase solid systems and materials such 
wood, paper, and clay low moisture 
quently, the diffusion equation satisfies drying rate data 
even though diffusion not taking place 
absence more detailed data, equations (10) and (11) are 
sometimes applied the second falling rate period 
However, neither type equation can used for 
drying rate curve like N-2 Figure 


EXPERIMENTAL 

Tray dryer 

The tray dryer was constructed some years previously 
for routine drying various materials and throughout 
these experiments was being used dry small experi- 
mental lots aluminum soaps. Air was circulated with 
in. vaneaxial fan directly connected motor 
mounted outside the dryer (Figure 2). one series 
measurements, air velocities between the trays were 
determined with Alnor Thermoanemometer and, over 
the height the dryer, were found vary parti- 
cular pattern from 4.8 6.5 ft./sec. with average 
Air was exhausted from the dryer with 
small blower and make-up air allowed enter, shown 
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Figure 2—Drawing the tray dryer used. 


Figure dry bulb temperature the air was 
automatically controlled regulating the flow steam 
heater and the wet bulb temperature regulating 
steam spray. 


Inside the dryer, the two tiers each had shelves 
spaced apart and each shelf held two in. in. 
trays. The upstream tier and the back half the down 
stream tier were not used but were filled with loaded 
trays that the air flow would resemble more closely 
the flow actual dryer. For the same reason, only 
alternate trays the remaining portion were used. Four 
trays aluminum soap were usually being dried and the 
tenth position from the bottom was used for the test 
tray. The regular drying trays had wire mesh bottoms 
but the test tray and the tray above had solid 
bottoms. The test tray was made stainless steel and 
was in. was fitted with twelve thermocouples 
placed line down the middle the tray parallel 
the air flow. Two couples soldered the bottom 
the tray and two placed and in. above 
the bottom. Wires couples were brought 
from below the tray through two rows small tubes 
in. apart and stretched across the intervening space. 
The tops the tubes were sealed with Khotinsky 
cement. The couples were connected two point 
Jones connectors, the otlier halves which 
tached the framework supporting the tray. Loops 
flexible copper wire connected the plugs thermocouple 
extension wires. The plugs and the soldered connections 
were the same level and only few inches apart 
that they would the same temperature. Temper- 
atures w ere recorded Ww ith a 12 point potentiometer re- 
corder. 


Two vertical rods extended from the tray supporting 
framework through the roof the dryer and the floor 
above it. horizontal bar attached the beam 
balance was clamped these vertical rods. 


Automatic weight recorder 


balance was modified considerably and pneumatic 
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Figure 3—Schematic drawing the weight recorder. 


components attached that the weight the tray 
could automatically recorded. The system used was 
essentially the same controller. Referring 
Figure decrease weight the tray moved the 
baffle away from the nozzle, causing the 
pressure decrease. This increased the relay valve out- 
put pressure hich restored the beam its original 
position. The relay valve output pressure was also con- 
was reversed that scale reading varied directly with 
weight. The range the instrument could changed 
moving the bellows unit along the beam and the zero 
set changing the tension the bellows spring, 
(2) altering counterweights either end the beam, 
(3) adjusting the nozzle position, (4) for 
slight adjustments, changing the position the recording 
pen. The instrument was set that zero reading was 
obtained with the empty tray position and was 
calibrated placing known eights the tray. During 
calibration, readings gave standard deviation 0.035 Ib. 
(for instrument range about but actual 
use, recorder weights the end run were error 
about 0.25 Ib. The recorder eights were always high, 
possibly because frictional resistance. 


Materials 


The materials dried were two types: aluminum 
soaps prepared here and commercial diatomaceous silica. 
One type aluminum distearate different 
aluminum dilaurates were prepared combining the re- 
spective sodium soaps with aluminum sulfate and sodium 
Freshly prepared material was used for each 
Physical properties the diatomaceous silicas are 
given Table Particle size data were taken from bulle- 
tin published anadian Johns-Manville Company from 
whom the materials were purchased while the bulk density 
and moisture content were determined here. These solids 
were mixed with water double conical solids blender 
and used several times before discarding. 


Procedure 

The test tray was placed the dryer, the thermo- 
couples plugged and the weight recorder set read 
zero. known weight wet solid was then placed 
the tray and spread evenly possible. The fol- 
lowing morning, later the day necessary, the tray 
was removed and the contents eighed. Samples the 
dry solid were taken for determination moisture con- 
tent with infrared moisture balance using 125 watt 
infrared lamp mounted 1.5 in. from the sample. 


(222 | | 
\ | 
FAN | 
STEAM 


TABLE 


Particle size, Hyflo 
microns Filter-Cel Super-Cel _Celite 545 

10-6 19.0 4.5 

Bulk density, g./ml. 0.22 0.28 0.33 


Aluminum stearate was dried three different temper- 
atures and the other materials two. Two runs were 
made with distilled water the tray instead wet 
solid. 


RESULTS 


Data processing 

Hourly values the average temperature the two 
top thermocouples, which were always well above the 
surface, were averaged give the air temperature for 
each run. The difference between these two couples was 
usually the order 2°F. Readings during typical run 
are plotted against moisture content Figure The 
weight dry solid, was determined from the weight 
and moisture content the end each run. The initial 
moisture content was computed from this and the original 
weight. Tray weights various intervals were read from 
the chart and the initial and final readings compared with 
actual any error was assumed have been 
linear function ing time. Corrected weights were 
then used calculate incremental drying rates which 
were plotted against time. Smoothed values were read off 
and plotted against moisture content shown Figure 
for one run. Values the rate fixed moisture con- 
tents were taken from this plot and tabulated. These 
rates for number runs with the same solid the same 
temperature were average and correlated with moisture 
content described later. 


Early runs 


Before correlating method had been derived, num- 
ber additional measurements and calculations were 
made. The diffusivity, was calculated from equation 


or TEMPERATURE 


w 
16410 
« = 
| 
= 
2 
~ € 
| 
2.10° 
e % > 
° > 
2 
wee = } 2 
° o* 6 
06410 
gé 
ay wee 
x > 
se 
> 
4 
2 6 2 


MOISTURE CONTENT, ory sone 


Figure for typical run with aluminum stearate. 


(9) and plotted against shown Figure for one 
run. For this plot, was assumed constant the 
final value. Assuming linear function the 
free water content gave different, but just 
content but plots for several other runs gave variously 
shaped curves correlation between and was 
possible. Thus, the simplified diffusion equation could not 
applied the data. 


The mean temperature the bed, which also shown 
Figure was arrived averaging the readings 
the thermocouples covered the bed. For this particular 
run, the bed depth was 0.68 in. the start and 0.30 in, 
free water content, was estimated that the thermo- 
couples half inch above the bottom became uncovered 
moisture content about 120%. Thus, six readings 
were averaged moisture contents above this value and 
four below it. plot the temperature any one point 
the bed would resemble the mean temperature plot 
shape. expected, the temperatures the bottom were 
generally higher than points near the middle the bed. 

The overall heat transfer coefficient was calculated 
from 


(12) 


at 
where 0.25 the estimated heat capacity the dry solid. 
The rate temperature change, was obtained 
from smoothed plot against all runs 
for which was calculated, showed steady 


Evaporation water 


Data for the two runs which water was evaporated 
were follows: 


Run Run 
Air velocity, ft./sec. 5.25 5.42 
Air temperature, 131.9 132.1 
Air humidity, 0.021 0.021 
Specific heat, 0.249 0.249 
Constant rate, 0.252 0.252 
Bottom temperature, 97.5 96.6 
Water temperature, 97.5 96.6 
Area bottom and sides, 4.70 4.70 
Area evaporating surface, 3.50 3.50 


The temperature the bottom the tray above the 
test tray was taken the same the air temperaure 
and the radiant heat transfer coefficients were calculated 
from equation (5) using emissivity and 
surface temperature 96°F. The convection 
heat transfer coefficients the the test tray 
were calculated from equation (2). 


Run Run 

Radiant heat transfer, 1.26 1.26 
167.2 167.8 

Bottom heat transfer, 1.52 1.52 
246.0 254.0 
Total heat transfer, 916.8 916.8 


These data were substituted into equations (1) and (8), 
which were solved simultaneously with 
humidity curve obtain the following: 


Run Run 

Convection heat transfer, 3.95 3.91 
503.6 495.0 
Mass transfer coefficient, 15.8 
Surface temperature, 95.1 94.0 
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The convection heat transfer coefficient higher than 
the 3.0 predicted equation (4). Perhaps there was addi- 
air turbulence because the water level was below 
the top the tray. 


Aluminum stearate 
Average drying rates aluminum stearate selected 
moisture contents are given Table for one temper- 


TABLE 


RATE ALUMINUM STEARATE 


Moisture, 

Number Rate, Standard rate, Error, 
solid (sq. ft.) (sq. ft.) 


0.0246 0.0055 


0.0717 0.0120 


ature. The air temperature ranged from 131.5 135.3°F. 
with average and standard deviation 
0.9. There was total runs but they did not all 
extend the higher moisture contents fewer points 
were available moistures above 50%. Drying rate 
plotted against moisture content for this and two other 
temperatures the lower portion Figure The 
straight lines the upper portion were obtained 
subtracting constant from the moisture content. This 
constant, which was found trial and error 0.01, 
represents the equilibrium moisture content. The straight 
lines, therefore, are plots drying rate against free 
moisture content, and the equations are 
the form 


FREE MOISTURE CONTENT, Ib ib. dry solid 


133.7°F 


1b. /(he)(sq.ft) 


ORYING RATE, 


ORYING RATE, 


MOISTURE CONTENT, dry solid 


Figure rate aluminum stearate. 
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Figure 6—Linear plot drying rate aluminum stearate. 


Figure arithmetic plot the same data and 
equations showing that the lines are concave downward. 
Values the constants and determined the method 
least squares are given Table can seen 


TABLE 


RATE EQUATION CONSTANTS FOR 


Number Equilibrium 
moisture, 


0.01 0.1989 


Figures and the equations fit the data very well. 
Correlation coefficients were above 0.998 and the standard 
errors estimate were 0.8, 7.6 and 3.7% respectively. 
These errors indicate the accuracy with which the regres- 
sion equations fit the average rate data while the standard 
deviations Table indicate the uniformity the rates 
from run run. The calculated rates 
errors Table will explained later. 

expected, straight line was obtained plotting 
log against the reciprocal absolute temperature but 
although decreased with temperature simple straight 
line relationship was apparent. 


Other materials 


Rate equation constants for six additional materials 
each two temperatures are given Table The 
effect temperature was reasonably consistent but 
temperature for two materials and 
decreased for the remainder. 


General correlation 


Since all values were between 0.54 and 0.72 was 
thought that constant value could used. This would 
make the equations less accurate for any particular solid 
but they would more generally applicable and only 
one constant would required for each combination 
material and temperature. additional advantage 
using constant was found that the resulting values 
gave more consistent slopes when log was plotted 
against the reciprocal absolute temperature. The aver- 
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TABLE 
RATE EQUATION CONSTANTS 


Air Number 
temp., 
runs 


Material 


Aluminum 
laurate 


laurate 


Aluminum 
laurate 


2330 


Filter-Cel 


Hyflo Super-Cel 1763 
2774 
3168 


Celite 545 


this was accepted the exponent for all rate drying 
equations. Thus, drying rates were given 


For each solid and temperature, was calculated 
various moisture contents using experimentally determined 
rates like those given Table for aluminum stearate 
133.7 Calculated drying rates each moisture content 
were then determined substituting the average rate 
constant into equation (15) and these rates compared 
with the experimental rates. Average percentage errors 
and average values are given Table which 
summary experimental conditions and results. 

The average drying rates Table were calculated 
from integrated form equation (14) and are the 
average rates drying from free moisture content 
150% the equilibrium moisture content. Drying rate, 
and moisture content are related 


TABLE 


CONSTANTS FOR TEMPERATURE CORRELATING EQUATIONS 


Average 
slope 


Individual 
constants 


Aluminum stearate 


Aluminum laurate —4,043. 
Aluminum laurate 
Aluminum laurate 
Filter-Cel 

Hyflo Super-Cel 


Celite 545 


—4,260. 
—6,526. 


and when this substituted into equation 
resulting equation can integrated give 


= — Xc -b 


where the time required reduce the moisture from 
The amount water evaporated during this 
time 


0.4831 


TABLE 


SUMMARY EXPERIMENTAL RESULTS 


Bed depth, in. 

Air Number 
temp., 
runs 


Material 
Initial Final 
Aluminum 144 24 0 0.36 

166 10 0 (). 28 


J\luminum 144 0 0. 
laurate A 17 4 0 0: 


Sluminum 4 0 62 0 


Filter-Cel 


0 0 
Supe! Cel 0.5% 0.4 


loading, 
ft. 


0 
0 
0 


0.7 


0 


0 


0.5 


0 


0.3 


0 
0 


Average rate 
moisture, Average 1.50 
dry solid (sq. ft.) 
0.0540 
1603 0.0774 
2102 0.1015 


Equilibrium 


0.01 
0.01 
0.01 


0.0632 
0.0944 


0.01 
0.01 


1853 0.0895 


0.01 
0.01 


1745 0.0843 
2908 0.1405 


0.015 
0.015 


0.015 1435 0.0693 
0.009 


0852 
1233 


0.002 
0.002 


0943 


0.004 


RATE CONSTANT, o 


RATE CONSTANT, o 


Material 
0.6253 1,190 
0.6691 
1,444 
0.6022 11,650 1,760 
Fig 
the 
onl 
anc 
Ta 
ene 
Fig 
69 
60 
49 62 
The Canadian Journal Chemical Engineering, April, 1959 


— 
z 
< 
CELITE 
z 
° 
HYFLO SUPER-CEL 
ALUMINUM 
LAURATE C 
« 
ALUMINUM 
LAURATE FILTER-CEL 
ALUMINUM 
LAURATE A 
STEARATE 


170 160 150 130 
TEMPERATURE, 


Figure 7—Variation rate constant with temperature. 


16) Effect temperature 


The rate constants Table are plotted against 

obtained with aluminum stearate was assumed hold for 

the other materials that were dried two temperatures 


17) only. The equations the lines are 

. . 
his and the constants and for each material are given 


Table The average slope corresponds activation 
energy 10,680 mole 5932 cal./g. mole 


water. This average slope and averages the respective 
pairs and were substituted into 
obtain more uniform set values (Table 6). 
Plots this equation for each material are given 
Figure Graphically, these lines have slope equal the 
Ll 
020 
3 ALUMINUM 
LAURATE 
HYFLO 
a SUPER -CEL 
FILTER-CEL 
ALUMINUM 
LAURATE A 
ALUMINUM 
STEARATE 
170 160 150 140 130 
TEMPERATURE, °F. 
Figure 8—Variation rate constant with temperature 
using average slope the lines Figure 
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average slope and pass through the mid points the 
respective lines Figure 


For the three diatomaceous silica materials, log gave 
almost straight line when plotted against bulk density 
but aluminum soap results were some distance from this 
line and not consistent although the general trend was 
the same. Log for the diatomaceous silicas gave 
smooth curve when plotted against average particle size, 
but, unfortunately, particle sizes the aluminum soaps 
were not available. 


Effect tray loading 


investigate the effect tray loading, rate constants 
were calculated for each aluminum stearate run and 
plotted against tray loading, final bed depth, and bulk 
density. the highest temperature, there appeared 
significant correlation but the two lower temperatures 
was obvious that there was none. The same was true 
when both slope and intercept were calculated for each 
run. general, correlated considerably better than 
giving significant correlation two temperatures but 
not the third. the lowest temperature, there was 
correlation between and drying rates 0.05 
but there was 0.50; the highest temperature 
the reverse was true. Thus appears that tray loading 
had significant effect over the range used, which was 
from 0.35 1.15 The drying rates the 
diatomaceous silicas substantiated this conclusion there 
was difference between rates tray loadings 0.4 


and 0.8 
CONCLUSIONS 


All experimentally determined drying rates were 
the falling rate zone constant rate period was 
observed. Moreover, since the overall heat transfer co- 
efficient decreased throughout each run, they were all 
the second falling rate zone, which characterized 
evaporation below the surface. Plots drying rate 
against moisture content were concave downward agree- 
ing with the conclusion Pearse, Oliver and Newitt 
that such curves are obtained with fine powders. 
number runs, there was tendency for the curve 
change concave upward the higher moisture contents. 
According Oliver and Newitt this probably happened 
because some water was being transferred capillary 
forces. 

experimentally determined drying rates were 
satisfied empirically derived equation that 
the commonly used linear equation (10) except that 
the moisture content exponent was 0.625 instead unity. 


Nomenclature 


Area evaporating surface, sq. 

Aw = Area of bottom and sides of the tray, sq. Ut. 

Constant (equation 14) 

Constant (equation 14) 

Constant (equation 21) 

Heat capacity humid air, dry air) 

Constant (equation 21) 

2.7183 

Mass velocity, (hr.) (sq. 

Saturation humidity the temperature the surface, 
dry air 

Convection heat from air the 
surface, (sq. ft.) 

Convection heat transfer air bottom and 
sides the tray, (sq. ft.) 

Radiation heat transfer coefficient from surtace above the 
the solid, (sq. ft.) 

Heat transfer factor detined equation 


Constant (equations and 11) 

Thermal conductivity the air, B.t.u./(hr.) (sq. ft.) 

Mass transfer coefficient, (sq. ft.) (unit 

Drying rate during the constant rate period, 
{sq it.) 

= Jotal heat transferred to the tray, B.t.u./hr. 

Heat transferred the solid, convection, B.t.u./hr. 

Heat transferred through the bottom and sides, B.t.u./hr. 
Heat transferred radiation, B.t.u. ‘hr. 

Reynolds number, ZVp dimensionless 

Temperature the solid surface, °R. 

Temperature the air, 

leniperature of the bottom of the tray, 

= Mean temperature ot the solid, °F 

Temperature the radiating surface, °F. 

Overall heat transfer coefficient between the air and the 
solid, B.t.u./(hr.) (sq. ft.) (°F.) 

= Air velocity, ft./hr 

Weight water the tray, lb. 

Weight dry solid the tray, 

Moisture content, dry solid 

Critical moisture content, dry solid 

Free moisture content, X*, dry solid 

Equilibrium moisture content, lb./lb. dry solid 

Emissivity the drying surface, dimensionless 

= Latent heat of vaporization, B.t.u./lb. 

Air viscosity, lb./(hr.) (ft.) 

3.14159 

Air density, ft. 


lime, hr 


czy 


to 


16) 
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Use Lead Chemical 


and Metallurgical Plant 


and 


WILLIAMS: 


The use lead the chemical and metallurgical 
operations Trail, B.C. described. 

wide variety corrosive environments are 
encountered, including phosphoric 
acids, sulphur dioxide, fluorides, and ammonium 
sulphate. found that except the presence 
fluorine compounds, antimonial lead more 
resistant corrosion than soft lead. 

The mechanical properties lead are reviewed 
with particular reference its high gravity, 
low mechanical strength, creep and fatigue proper- 
ties. Principles plant design give adequate sup- 
port, provide for thermal expansion and contrac- 
tion, and minimize the effects creep and fatigue 
are enunciated, and illustrated proven designs 
that have given long periods satisfactory service. 


has history usefulness dating back anti- 
quity. our present industrial economy, flooded 
with new custom-built products, each designed its 
specific application, lead has survived because its unique 
properties make still the first choice, and often the only 
choice, for many services. Its use chemic: 
plant construction dates back the very foundations 
chemical industry, and reported have been used 
phuric acid chambers. 

the subsequent development the Chamber Pro- 
cess for sulphuric acid manufacture, has remained the 
only material construction with the 
sistance corrosion which can easily fabricated into 
such large units. 

merits have gained general recognition, and 
experience the use lead chambers and their auxiliary 
equipment has provided foundation for much wider 

Outstanding, but probably less obvious advantages 
construction its versatility and its ease 


IManuseript received July 31, 1958 

2Research Engineers, Research and Development Division, The Consoli 
dated Mining and Smelting Co. of Canada Limited, Trail, B.C 
’Superintendent of Development, Engineering Division, The Consolidated 
Mining and Smelting Co. of Canada Limited, Trail, B.C 

Contribution from The Consolidated Mining and Smelting Co. of Canada 
Limited, Trail, B.C. and presented at the pit AICHE-CIC. Chemical 
Engineering Conference, Montreal, April 20-23, 1958 
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forms including large sheets, large 
and variety cast shapes. easily cast the 
simplest foundry techniques. Pipe and sheet are ‘easily 
and welded, and they require expensive plant 
workmanlike job. Piping and linings can 
cut, patched altered the use simple tools and 
with little difficulty adapt plant changing require- 
ments. Salvaged plant can transplanted new loca- 
tions, where good new construction. Shop 
scrap can easily remelted and cast pipe 


hangers fittings simple chill 


addition these merits, lead equipment requires 
upkeep painting because the bare metal remains neat 
and smart. When finally scrapped, the scrap value 
the lead equal substantial proportion the 
first cost. 


would seem inevitable that any material construc- 
tion with many would sutter some 
handicap. the case lead, the and ductility 
inherent mechanical weakness, which has been 
aging its wider use. reviewing the cor- 
rosion-resisting and the mechanical properties lead, 


have been developed to overcoihe this weakness. 


Lead corrosion resisting structural material 

Comparatively few types alloys lead are used 
for chemical There are only two 
major ones. “Soft” lead may relatively pure lead, 
containing Cu, known “Chemical 
type “hard” lead, normally lead-antimony 
the range antimonial lead being the most 
common. 

required, pure lead (99.994 Pb) generally used. 
lead equipment, and especially antimonial lead, for use 
within its plants shown the statistical data 

Six percent antimonial lead used preference 
soft lead for two reasons: (1) because has substantially 


ed, 

Book 
976 

1096 

936). 
805 
them, 


LEAD EQUIPMENT FABRICATED ANNUALLY 
(5-Year Average; Tons) 


Form Lead** 
_Lead 

Pipe 160 
Castings Nil 110 


All torms 100 


*Exclusive of silver-lead anodes and associated equipment for the electrolytic 
zinc process, equivalent to 725 cons/year 

**Refers to both pure lead (99.99+% Pb) and “Chemical Lead” (Pb + 
006% Cu) 

** Large proportion of soft lead sheet used for flashing 


better physical and (2) because has 
equally good better corrosion resistance. The mech- 
anical properties lead and antimonial lead alloys will 
described later. That the corrosion resistance anti- 
monial lead generally equal better than lead 
somewhat contrary popular opinion. However, 
will shown the followi ing discussion, this the 
general conclusion based corrosion testing work and 
long practical experience the field. This observation 
course, restricted the corrosive media encountered 
Trail, but the sulphur dioxide recovery system, sul- 
phuric acid plants (over total Canadian produc- 
tion) wet phosphoric acid plants, and 
electrolytic zinc (zinc sulphate) process cover wide 
range conditions. 

The generally good corrosion resistance the 
antimonial lead alloy may best illustrated consider- 
ing specific plant corrosion tests. Trail, Cominco 
operates phosphate fertilizer plant based the Dorr 
wet process for the production phosphoric acid and 
lizers. This plant and process have been described 
the literature 

With its various combinations phosphoric, sulphuric 
and hydrofluoric acid solutions and slurries, all 
latively high (60-120°C.) operating 


TABLE 


CORROSION RATE VARIOUS METALS EXPOSED 
ACID SLURRY 
Spool exposed for days 


Material Corrosion Rate (ipy) 

Lead. . 0.0054 
AISI 316L. 0.0062 
0.0060 


presents most difficult corrosion problem. 
due the omnipresence the sulphate ion, lead 
useful material construction for many parts the 
plant. Six percent antimonial lead used almost ex- 


Fortunately, 


clusively throughout this plant. For slurries, 
the antimonial lead alloy somewhat better than 
pure lead, shown the data Table course, 
where abrasion serious problem, lead equipment must 


The results field tests carried out the various 
mixtures phosphoric and sulphuric acid resulting from 
the countercurrent reaction and filtration process em- 
ployed have shown that neither lead nor its alloys are 
corroded these solutions. Test specimens pure 
lead and antimonial lead exposed mixtures 
phosphoric and sulphuric acid ranging composition 
from H,PO,, and temperature from 
show negligible weight losses after days. 


Comparative tests lead, antimonial lead and 
AISI Type 316 stainless steel (Table carried out 
various sections the plant ventilating and fume re- 
covery systems serve illustrate the general conclusion 
that dilute phosphoric- sulphuric acid and ammonium 
phosphate solutions the corrosion rate antimonial 
lead the same order, lower than, pure lead. 
should noted here, however, that wherever fluorine 
compounds, (mainly and tend concen- 
trate, the corrosion rate lead, and for that matter most 
other metals, increases rapidly. 


This concentration fluorine 
the only exception the general conclusion that 


LEAD AND STEEL EXPOSED VARIOUS VENTILATING DUCTS AND SCRUBBERS, 
PHOSPHATE FERTILIZER 


(All corrosion rates 


Fertilizer Section Reaction Agitator System 


Material 
Gases 
Reaction Agitators 


Scrubber 


trace 70°C. 


50°C. 


Fertilizer Dust Acid Section 


Recovery Reaction 
Agitators, 
Exit Scrubber Ventilation 
From Scrubber 
Negligible 0.018 
0.003 0.013 
0.0001 Negligible 0.3 


(slight 


Spool exposed days. 


Spool exposed days. 


Note \ir and steam plus phosphoric acid, mono-ammonium phosphate, trace SiF 90°C 
Air and steam plus sulphuric acid, ammonium 
\ir and steam plus phosphoric acid, fluorine compounds, 


Spool exposed days. 
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TABLE 


Corrosion RATE VARIOUS METALS EXPOSED 
FLUORINE REMOVAL EQUIPMENT 
PHOSPHATE FERTILIZER PLANT 


aerated and agitated. 100°C. Spool exposed days 


Material Corrosion Rate (ipy) 

Hasteloy C........ 0.066 

0.12 
Carpenter 0.28 


Specimen Disappeared 
Chromel A.. Specimen Disappeared 
Hasteloy B... Specimen Disappeared 
AISI Specimen Disappeared 


Aluminum, 


antimonial lead has equal better corrosion resistance 
than lead this plant. 


Fluorine removed from the plant circuit mixing 
proportion the phosphoric acid produced with con- 
centrated sulphuric acid and agitating (stripping) the 
resultant hot solution with air. This operation, which 
undoubtedly the most corrosive single operation within 
the Trail plants, carried out especially built units 
called “strippers”. The strippers are made from cast 
antimonial lead sheet, one inch thick, carbon brick lined. 
With some difficulty, corrosion spool was exposed 
the conditions inside stripper with the results shown 
Table 


will noted that lead the only possible material 
construction for this service. actual fact, the cor- 
rosion rate both lead and antimonial lead, while satis- 
factorily low, may vary from place place inside the 
stripper, depending fluorine concentration and degree 
aeration and agitation. practice, the strippers are 
fabricated from antimonial lead order take 
advantage the slightly better (at the operating temper- 
ature) mechanical properties this alloy. 


The well-known corrosion resistance lead and its 
alloys gases and solutions containing sulphur dioxide 
utilized many ways for handling the numerous gases 
and solutions this type processed Trail. Tables 
and illustrate pical results. Again, the corrosion 
resistance antimonial lead these solutions has 
been found equal superior that pure 


lead. 


TABLE 


CoRROSION RATE VARIOUS METALS EXPOSED 


Saturated Water plus 45°C. 
Days’ Exposure 


Corrosion 


Material Rate (ipy) Remarks 


Sb-Lead Negligible visible attack 


Lead Negligible visible attack 
AISI 316 Slight general attack 
304 0.007 Slight pitting 

AISI 310 Pitted 

AISI 410 Specimen disappeared 
AISI 446 Specimen disappeared 


Mild Steel 
Aluminum, 


Specimen disappeared 
Specimen disappeared 


The Canadian Journal Chemical Engineering, April, 1959 


TABLE 


RATE VARIOUS METALS EXPOSED 
Jet CoNDENSER, AMMONIUM SULPHATE PLANT 


(0.1 gpl 0.01 gpl SO:, Turbulent Flow, 
Temperature 40°C. 
Spool Exposed for days 


Corrosion 
Material Rate (ipy) Remarks 
Sb- 0.003 Darkened 
0.001 Appearance unchanged 
AISI 347... 0.001 Appearance unchanged 


Mild Specimen disappeared 


Cast 0.24 Severe pitting 
Wrought Iron. 0.55 Pitted 
Aluminum, 0.008 Some pitting 
0.04 Pitted 

Brass 0.05 Slight pitting 


The data presented Table are particular interest. 
The Doyle Scrubber requires that the gas stream 
impinged high velocity into liquid. Normally the 
incoming gas stream carrying some solid burden direct- 
nozzle into water. The outside the nozzle 
exposed the water, saturated gas and splash. While 
lead not normally thought abrasion resistant 
material, practice, lead has proven 
admirably well suited for handling 
suspensions. 

Although the corrosion rate for Type 316 stainless 
steel was negligible the test results shown Table 
pitting type corrosion will occur those areas 
where sludge build-up experienced. Lead preferred 
material construction this environment because 
shows tendency corrode pitting. 


Historically, noted above, the first important use 
lead the chemical industry was the production 
sulphuric acid the ch: amber process. this process, 
the sulphuric acid contains some oxides nitrogen. 
the search for improved lead alloys for chamber plant 
construction when the process first came into large scale 
use the early was the practice the day 
test new alloys the laboratory hot strong acid. This 
technique culminated the so-called 


the flash test, specimen lead lead alloy 
specific dimensions was immersed 
phuric acid and the temperature the acid was raised 
prescribed rate. certain temperature, precipitate 
lead sulphate rapidly formed the acid. This tem- 
perature was defined the flash point. differ- 
ent impurities and alloys produced different results. The 


RATE VARIOUS METALS EXPOSED VENTILATING 
GASES FROM THE LEAD SMELTER SINTER CRUSHING PLANT 
Spool exposed for 140 days the gas stream inside Doyle 
Scrubber. Gases contain some solids, trace are saturated 
with water vapour, and range temperature from 40°F. (winter) 
(summer). 


Corrosion 


Material Rate (ipy) Remarks 
Lead 0.0005 Some slight mechanical damage 
316 Negligible apparent change 
Mild Steel 0.07 Severe attack 


apparent change 
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TABLE 
RATE VARIOUS METALS EXPOSED 
SULPHURIC ACID 
Acid Aerated and Agitated Temp. 25°C. 

Spool Exposed for days 


Corrosion 


Material Rate 


<0.001 


assumption was that the lead with the highest flash-point 
would give the best service. the course time, 
was observed that the “order-of-merit” indicated 
the flash test did not correlate with the results actual 
corrosion tests carried out the field 
curately simulated service conditions, and the test was 
discarded. Today, the possible short-comings 
quick laboratory test and the merits long term pro- 
perly supervised corrosion tests under service conditions 
are widely recognized. 

One the alloying additions that gave consistently 
low was antimony. Also, short term 
tests hot concentrated chamber acid, antimonial lead 
alloys gave somewhat higher weight losses than other 
and alloys. For this reason the antimonial lead 
alloys were first erroneously regarded unsuitable 
for other types chemical plant construction. For- 
tunately, actual tests have shown that the 
strength alloys have satisfactory cor- 
rosion resistance many media, and these useful alloys 
have found increasingly ider use. 


Another factor contributing the increased use 
antimonial lead alloys has been the increasing importance 
the modern contact sulphuric which 
brought about fundamental change the type equip- 
ment and materials required for acid plant construction. 
the five contact sulphuric acid plants Trail, iron 
and steel are commonly used for handling concentrated 
sulphuric acid, and brick-lined steel for handling hot 
concentrated acid. Both pure lead and antimonial lead 
are used for handling the more dilute sulphuric acid 
solutions. While many soft lead structures have been 
service for many years, modern practice 
antimonial lead for most construction. That 
monial lead may used for handling even concentrated 
sulphuric acid (cold) shown the data Table 


The mechanical properties lead 

differs from the common structural 
asmuch its mechanical properties ordinary temper- 
atures tend toward those plastic solid rather than 
elastic rigid one. For this reason, laboratory results 
the strength, creep, and fatigue not have 
the usual significance when interpreted for 
poses. Below given summary the main results 
studies these properties, and shown 
later that with proper design, results can obtained 
which are more satisfactory than would anticipated 
from the 

Widely varying values have been obtained for the 
tensile strength, yield point and hardness lead. These 
values are influenced not only minor variations 


esults laboratory tests. 


position, but also the actual conditions the test. 
Thus, has been that figure recorded 
any value unless the rate application the load 
also given. illustration this point, the following 
figures have been quoted for particular sample high 
purity lead different rates extension: 


Time fracture 1.5 minutes: ultimate tensile strength 

Time fracture minutes: ultimate tensile strength 

The tensile strength lead under short time test 
evidently little significance the design equipment 
for long life. More strength under 
sustained loads for long periods time, which de- 
signated the creep strength, and the strength under re- 
peated cycles stress strain, described the fatigue 
strength. 

Since the phenomena creep and fatigue any metal 
are more significant one approaches its melting point, 
lead has been the subject intensive research because 
exhibits properties ordinary temperatures which can 
only studied other structural metals very high 
temperatures. For details this work, reference ‘should 
made papers published Professor Moore 
and his associates the University Neill 
Greenwood Australia, Bailey and McKeown 
the British Non-Ferrous Metals Research Association, 
and others. 


Creep lead and lead alloys 


Moore, Betty and Dollins have quoted experimental 
evidence show that there limiting stress under 
either long-continued steady load, with cycles 
repeated stress, below which fracture lead will not 
eventually occur. has also been found that although 
with antimony increases the tensile strength 
pure lead about times, the increase creep 
strength only about three times ordinary temper- 
atures. Greenwood gives comparison 
this (Table 9). 


These figures suggest that alloy containing around 
antimony has long-life properties better than either 
pure lead higher antimony alloys, both ordinary 
somewhat higher than ordinary temperatures. Since 
antimony additions above make the lead difficult 
work, these data give for the general 
practice adopting the alloy standard material 
for lead construction. temperature raised, the 
improvement creep strength due the antimony de- 
creases and eventually disappears temperatures around 
120°C. (250°F.) See Table and Figure 


TABLE 


COMPARISON TENSILE STRENGTH AND CALCULATED BREAKING 
STRESS FOR 30-YEAR LIFE 


Strength 30-Year Strength 30-Year 
Life p.s.i. Life p.s.i. 
| | 
Sb 2,600 500 | 2,000 325 
%Sb 7,700 775 5,500 500 
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TABLE 
ALLOWABLE FIBRE STRESS EXTRUDED PIPE 


Maximum Allowable 


Temperature Equivalent Fibre Stress p.s.i. 
Lead Lead 
370 
104 180 340 
140 162 280 
158 153 254 
176 144 
194 136 195 
100 212 127 165 
110 230 118 137 
120 248 110 110 
130 266 100 
247 477 535 


where safe working pressure p.s.i. 

maximum allowable fibre stress from above table 

thickness pipe wall inches 

inside diameter pipe inches 
For mechanical structural reasons, sometimes wise 
allow additional wall thickness above that derived from this 
formula. 


attempt has been develop formula 
express the stress required produce failure lead 
alloys after the elapse any given period time, but 
the extrapolations have been questioned because the 
data did not fall straight line, and also because the 
creep pure lead has been shown discontinuous 
process deformation with alternating periods rapid 
and extremely slow extension 


Fatigue lead and lead alloys 

Studies fatigue lead and lead alloys have been 
largely related distinctive type 
cracking cable sheathing service. Similar failures 
have also been observed lead piping which has been 
subjected bending stresses immediately after fabrica- 
tion. The same type failure has been reproduced ex- 
perimentally use Haigh testing machine imposing 
repeated cles stress test specimen. Other types 
test have also been devised which impose repeated 
bending stresses the specimen, reproducing more close- 
the conditions encountered service. 


fatigue stresses decreased the temperature raised. 
Thus found that sample which failed 
million reversals 4°C. lasted only million 
32°C., while Beckinsale and 
similarly found that lead which failed 1.17 
versals Haigh machine ordinary temperatures, 
failed 0.988 million 100°C. and 0.556 million 
150°C. all tests being made under repeated stresses 
+560 p.s.i. 


The rate reversal cycles stress has 
portant influence the endurance under 
the number failure decreasing the frequen- 


fatigue stresses, 
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Allowable stress temperature curve 
for extruded 6% Antimonial 
lead 


Allowable stress 
for extruded soft 


ALLOWABLE STRESS 


100 200 400 500 
TEMPERATURE °F 


Figure fibre stress lead. 


decreases, though the actual time failure may in- 
crease (12, 13), 


The effect alloying elements generally cause 
grain refinement, hich has the effect increasing the 
fatigue resistance lead, but even alloys comparable 
grain size, has been found that the addition anti- 
mony, tin, other alloying elements increases the fatigue 
resistance. This improv ement persists the temperature 
raised (14, 15). 


From this brief review the creep and fatigue pro- 
perties lead, evident that even ordinary temper- 
atures, the mechanical properties metallic lead are 
difficult define. For plant design purposes, therefore, 
guidance must sought from empirical data based 
actual experience. 


series such figures developed for the design 
extruded pipe systems was published Roll and are 
given Table 10. 

Although laboratory tests show that any continuous 
loading lead, any cyclic stress will entually cause 
failure through creep fatigue, practical experience has 
shown that with proper design and adequate support, lead 
construction will last indefinitely. 


Lead tank and tower walls, for example, have been 
found undergo certain amount distortion 
bulging between their steel-frame support during their 
initial period service, and then remain unchanged for 
long periods time after arriving their stable form. 

Tanks are still service Trail which were installed 
years ago having had very little repair maintenance 
the interval. 


Design factors lead installations 


The successful use lead construction conditional 
the recognition the differences the properties 
lead which distinguish from the more common struc- 
tural materials. The most important these are, 
course, its mechanical weakness, its lack 
which results permanent deformation under stress, its 
tending creep under any sustained load, low 
fatigue limit under repeated eight adds 
first requirement any 
design is, therefore, provide adequate support. 


the problems design. 


moderately elevated temperatures, the problem 
support more exacting, the strength lead and 
its alloys falls off considerably with rising temperatures 
(see Figure 1). 


400 
| 
} 
| 
| 
calculate safe working pressure lead pipe, use formula 


second requirement design the provision for 
expansion and contraction with changing temperatures. 
The situation one which the installations are 
constructed such way that the lead free move 
without hindrance under the influence changing tem- 
peratures. free movement can provided then 
there stress that could result failure. the 
traditional lead chambers sheets were supported from 
their tops and hung free with lower edges over- lapping 
the up-turned edges the floor the chamber. These 
chambers were about high, which about the limit 
for free-hanging sheets lead. the requirements for 
support not permit free contraction and expansion, 
then the fixed points supports should not too wide- 
spaced, that the amount expansion taken 
each section remains Often, such design 
satisfactory and there tendency towards failure 
through fatigue creep. large installations, properly 
designed expansion joints are provided that excessive 
stresses cannot develop. 


The adoption antimonial lead with its greater 
strength and fatigue resistance has made the problem 
support somewhat easier solve, but has also introduced 
new problems. Being less ductile, less able 
adjust itself local stresses and somewhat more prone 
cracking points concentrated stresses. Once this 
recognized design can adjusted avoid such stress 
raisers and trouble encountered. 

The application these basic principles design 
lead construction, adequate support, proper pro- 
vision for thermal changes, and avoidance stress 
raisers illustrated number installations actual 
use the Trail Plants. Since the practice recent 
years has been toward the general use antimonial 
lead, the specific examples described have been construct- 
that alloy. The principles involved, however, have 
equal application soft lead construction and the 
designs can adapted that material recognition 
the differences its properties. 


Lead piping 

One the great advantages lead piping, addition 
its corrosion resistance, its adaptability 
ease with which can installed. When pipe bends 
are required, straight pipe may filled and packed with 
sand, and bent. During the operation, the outer wall 
becomes thinner and the radius the bend less 
than four times the inside diameter the this outer 
wall “heel” reinforced “burning on” added 
layer lead. Shorter radius bends may made 
means mitred joints between short sections pipe. 
Alternatively, variety pipe fittings may pur- 
chased stock items, cast chemical lead anti- 
monial These are supplied four ways 


With solid lead flanges welded place. 

With cast iron steel slip-on flanges. 

With over-length for flanging the field. 
With plain ends for welding. 


Hangers on 


alternate sides of 
hanger bor at 


Hard 


dia O 


Figure 2—Method hanging lead pipe. 


When flanges are used, the diameter and drilling nor- 
mally conforms American Standards 
quirements for Class 125 cast iron flanges. 


The traditional method supporting pipe runs 
continuous wooden supports steel troughs has the dis- 
advantage that any leakage causes damage the supports, 
that free expansion and contraction are hindered, and that 
the location any leak not always evident, easily 
accessible for repair. 


The preferred method hanging lead pipe 
Cominco plants illustrated (Figure 2). This method has 
the merit infinite adaptability, neatness, and negligible 
up-keep cost since lead does not require any protective 
painting. provides unhampered movement under stresses 
produced expansion and contraction, and this way 
the ideal support for lead. Since the supports are all 
above the pipe, they are immune from damage from the 
contents the pipe leakage should occur, and the 
point the leakage easily located and accessible for 
repair. 


The hanger bars are cast simple chill moulds the 
shop, using shop scrap. 


Pipe runs diameter can supported this 
manner. Lengths feet long have been prefabri- 
cated complete with hanger bars, hoisted position, and 
the joints welded into position with butt welds the 
field. Hanger bars are hard (6%) antimonial lead, 18” 
long, and are welded the pipe foot centres 
hard lead pipe. 


Lead sheeting 


Welded joints are standard construction chemical 
engineering practice where lead linings coverings must 
protect the supporting structure from infiltration cor- 
rosive gases solutions. general, lap joints (Figure 


For Strip asbestos 


Figure 3—Typical welded butt joint. 


are used for joints vertical sheet lead. horizontal 
surfaces, where the free movement the lead may 
restrained the weight the lead itself, the 
pressure liquid the lead, there tendency for 
the bend the lap joint act expansion bend, 
localizing the stresses and strains and causing failure along 
the lap the joint. For this reason, butt joints are used 
horizontal surfaces (Figure 4). 


NOTE When hydrostatic 
other pressure the 
weight of Lead itself ona 
floor acts ona 
Lead surface with lap joints 
without allowance for 
expansion or contraction, 
any expansion or contraction 
the Lead sheet will 
cause the bend at the lap 
joint to act as an expansion 
bend, localizing the stress 
and strain. Experience 

indicates that the Lead will eventually fail 

along Line"A”. In this application, butt 
joints are recommended with full provision 
for expansion. 


TYPICAL LAP JOINT 


lap joint. 
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Tanks 


Unlined wood stave tanks are often used for holding 
liquids without ill effects except for the corrosion the 
stecl hoops. Lead pipe has been used enclose the hoops, 
indicated (Figure 5). This system unsatisfactory 
because any injury the pipe allows entry the 
solution, will lie the lower half, and may corrode the 
stecl beyond repair before the damage apparent. Effec- 
tive protection provided draping lead strips over the 
pipe, leaving the hoop open for inspection, and allowing 
any drain away without causing damage the 
hoop. 


Since unlined wood stave tanks depend the swelling 
the staves for solution tightness there tendency for 
the staves dry and shrink when the tanks are lead lined. 
The tightening the hoops maintain structural rigidity 
would have most destructive effect the lead lining. 
avoid this condition, reinforcing ring wood 
constructed around the inside top the tank maintain 
its original diameter and shape. 


Details the design such tank are shown (Figure 

the height the wall does not exceed feet, 
support the vertical sheets not necessary. 
will also noted that the side sheets are continued 
beyond the angle ‘of the base and butt-jointed the floor 
lining circle about inches inside the circumference 
the base. This provides some freedom for expansion 
and contraction the wall lining without any concentra- 
tion stresses about fixed points, along fixed lines. 


installations over feet high, vertical seams are 
supported for their whole length. making the lap joint, 
the vertical edge the sheet first secured the wood 
with flat head nails about 10” apart, and clear 
the edge the sheet. The second sheet lapped well 
over the nail heads and welded give lap not more 
than 2”. The whole seam then protected and supported 
extruded lead supporting bar (Figure 
63” 3”) which secured carriage bolts intervals 
The heads the bolts are protected welded 
cover lead sheet. Earlier methods support depended 
the use steel strap through which the supporting 
bolts were passed, the being protected welded- 
cover sheet lead. Such method produced concen- 
trations stresses the sheets the vicinity the 
supporting bolts. The adoption the extruded support- 
ing bar, with continuous welded joint the sheets 
both sides, has eliminated this source weakness. 


Walls all tanks, where possible, are lined with con- 
tinuous sheets from top bottom. This avoids concen- 
trations stresses where over-hand burning lap joints 
crosses the vertical seam. sheets sufficient length are 
not obtainable from the supplier, two shorter sheets can 
butt-welded end-to-end, before erection. 


Flues 

the design flues, the problems support and 
expansion both become more difficult the cross-section 
and the the flues increase. 


Circular flues 

Circular flues diameter 42” have been 
constructed lead and are operating satisfactorily. 
the construction these flues, antimonial sheet lead 
about pound weight built around cylindrical form 
and provided with welded-on circular lead stiffener 
12” intervals. 
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Figure 7—Lead supporting bar. For use fabrication 
lead vessels stiffeners. 


All joints these flues are butt joints. joining 
sections larger flues, the lower half can welded from 
the inside and the upper half from the outside. 

flue supported lugs the stiffener rings, 
resting longitudinal beams either side the flue, 
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Honger spacing as - 
determined Dy design 
Single hanger for 
flues up to about 
16"diam 
Double hangers 
for flues over 
6"diam 
leod 
Reinforcing ring spocing 
as determined from design 
( usually about 12") 


CIRCULAR FLUES AND HANGERS 


Mord lead st tfeners 


Anchorlfbo 1s a? 
about 114-Oc-c 


RT 


3" Lead pipe split longitudinally 
The remaining half is again 
Split ond ali are welded to 

form section shown 
Bend to fit flue and 
butt weld 


Sections of flue made in 
7—6° lengths and butt 
welded together in fieid. 


Figure 8—CIRCULAR FLUES: 
(a) flues and hangers. 
(b) flue showing method support. 


Suggested method allowing for longitudinal expan- 
sion. 


\ Temporary 
staging 


weld 


= 
“12 Ib.Lead 


SECTION 


Expansion vends 


spaced ot GFI.c , L@p weld where necessary or 


use expansion bends at 

necessary points conforming 

== to available sheet 
| widths 


Support strips to be continuous 4 SECTION B—@ 
except at expansion bends 

Avcid wreaking weld ta main lead sheet 


& weid will Guilast an intermittent one 


flue with structural steel basket 
construction, 


hangers hooked into holes drilled the stiffeners, 
the former case, the lugs are welded bearing plate 
and anchored with bolts 4-foot intervals. This 
bearing plate absorbs some expansion stresses and pro- 
tects the flue from deformation. (See Figure 8). 


When supported hangers, single row, hooked into 
holes the top the adequate for flues 
16” diameter. For larger flues, double hooked 
into holes drilled the sides the stiffener rings are 
required. Hangars are provided feet apart give 
adequate support. (See Figure 8). 

flues which are provided with bends, often not 
necessary install expansion joints. But long straight 
runs where flue joins another structure, expansion 
joints soft lead, the type illustrated, are 
(See Figure 8). 


Towers and flues 


the past, was Cominco practice build towers 
and flues laminated wood construction lined with lead. 
These had the disadvantage that the sheets were not 
accessible for inspection from the outside, and they have 
recently been superseded steel basket type con- 
struction support the lead sheet. These have the 
advantage that they can easily inspected and repaired, 
sometimes while the vessel still service. 

Structural details steel basket construction are 
illustrated Figure Continuous pound lead strips are 
welded the sheet before erected, and these are 
anchored the horizontal angle sections the steel 
basket. most important that these strips continuous, 
otherwise point weakness developed the gap. 


Spacing to conform to pe ag leas NOTE -This construction was used 
nor when was foreseen that 
another similar flue may 
be required along-side. 
parts were lead 
covered to allow this wail 
5" Pipe columns to become a dividing wall 
covered with in the future. Variations 
10 Ib. Lead sheet service requirements 


PLAN SECTION A—A are conceivable. 


Expansion 
12 |b. Lead sheet 


Pipe at 
with lead pipe 


staging 


lead flue with pipe basket construction. 
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supports are continuous excepting where 
expansion bend inserted. 


later development this type flue construction 
the pipe- illustrated Figure 10. 
Expansion bends and vertical joints this type con- 
struction are combined the design and located the 
circumference the vertical pipe column. 


Lead clad steel 

plants reach their extremes large flue 
exposed differentials, vibration, fluctuating 
temperatures, corrosive gases condensa- 
tion moisture. All these conditions are encountered 
gas-handling system conducting hot, dust-laden gases 
from the lead sinter plant Cottrell Plant and sulphur 
dioxide absorption units Trail. 


diameter flue was constructed handle gases carrying 
heavy dust burden and containing about SO,, enter- 
ing temperature high 400°F. the discharge 
end the flue, booster fan forces the gases down 
cooling tower ft. high and ft. diameter, con- 
current with fine spray water which cools the gas 
and wets the dust particles condition them for treat- 
ment Cottrell precipitator. Corrosive conditions are 
increased when there condensation moisture during 
shutdown periods when the flue cools down. Stainless 
steel could not used because pitting due contact 
corrosion aggravated the presence fluorine com- 
pounds. 


meet these conditions, the whole the flue was 
constructed steel with lead-clad interior lining, and 
the cooling tower was constructed lead-clad steel, 
supported concrete sill and base. The steel casing 
the booster fan, received from the manufacturers, was 
lead spraved Trail, using Metco wire gun, and the 
spraved coating was further sealed with coating 
epoxy -phenol resin. The steel fan runner was not coated, 
since easily replaced and not subject condensa- 
tion moisture. 


The lead cladding was applied cylindrical sections 
the flue steel, units 9-ft. diameter and 8-ft. length. 
Internal surfaces these cylinders were cleaned, 
and then given first coating lead rotating 
shallow bath, that the surface the steel dipped into 
the molren Jead. lead adhered the unfluxed external 
steel surface. The main coating lead, about ft., 
(5/32 in.), was then applied using equipment and methods 
developed Knapp Mills Incorporated Long Island 
New York. The sections were then welded into 
units about and the internal lead cladding 
applied over the welded joints hand burning. These 
sections were provided with end flanges and were trans 
ported the site and erected. The flanges were bolted 
up, and the cladding over the joints “burned” 
In position. 


The cooling tower was constructed lead clad mild 
clad the Knapp Method, supported concrete 
construction, each plate being spot welded form 
cvlinder, then cleaned, fluxed and clad the same way 
the sections the flue. plates were then re-rolled 
the correct curvature and welded ‘The 
cladding the welded joints was then sealed cladding 


lead over them hand methods. vertical walls 


the Jower were protected Ib. (3/16 


in.) sheet lead welded vertical supporting strips lead 
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inserted the concrete. The floor was protected sheet 
lead laid over the concrete, after the application 
intermediate dressing asphalt. This lead was protected 
course acid proof brick, against the possibility 
mechanical damage cleaning out accumulated sludge 
should necessary. 


Finally, the external steel surfaces the flue and the 
cooling tower were sandblasted and primed with zinc- 
rich paint provide cathodic protection, then given 
outer coat silicone-base heat resistant paint. 

Lead cladding homogeneous bonded lead prov ides 
the ultimate support lead construction. Steel the 
normal base material but others can used. Operating 
temperatures can carried close the melting point 
lead, thus extending lead’s usefulness into wide new 
field. Hand burning, pouring (casting) the machine 
methods Knapp Mills may employed. 
Rotating elements such fan runners, high temperature 
steam heating coils, and lined pipe for high pressures are 
easily constructed. 


Due the difference the temperature coefficient 
expansion between lead and steel, there long term 
tendency toward failure cracking when 
extreme and frequent temperature fluctuations are en- 
countered. case penetration the lead any one 
spot, there creep the corrosive media behind the 
lining and repairs are easily effected. Good workmanship 
essential, and the few reported are 
ascribed discontinuities lead and lack com- 
plete bonding. 


Economics 

illustrate cost lead work, Table gives con- 
struction estimates typical tank several corrosion 
resistant materials. tank ft. diameter and ft. 
high for storage acid solution, has been used for the 
example. Labor was figured shop rate 
$4.00/hr, sheet lead and extruded lead the current 
quoted price and other materials current 
prices. 


botal 
Cost 


lead lining wood staves 


rubber lining '4” mild steel 


lead cladding mild steel 


stainless steel (average plus) 9000 


design each case normal tor the service, and 
the intention merely show the relative costs 
construction lead and alternative materials. 


competitive with other materials 
serap value the lead quoted, would 
represent an additional $000 (at current serap quotations 
which would recoverable the end 
the life the installation. complex 
the versatility lead sheet construction would possibly 
rather lower costs relation other 
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During the service life the installation, the low 


maintenance costs, the excellent corrosion resistance, the 
ease and reliability with which modifications can made 
lead sheet equipment, and the excellent service 
high temperatures given lead-clad plant, contribute 
further the overall economic advantages lead con- 
struction. 
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The Pyrolysis Heavy Crude Oil 
Fluidized 


and 


The pyrolysis heavy high sulphur crude oil 
inert fluidized bed temperatures between 
550 and 800°C. has resulted the production 
olefin rich gas and product containing 
increasing amounts aromatic material 
creasing temperature. 


Analysis the process gas has shown that for 
average reactor residence times more than one 
second, equilibrium approximated for the system 
ethane-ethylene-hydrogen. There 
the data for the analogous system propane- 
propylene-hydrogen but for most conditions signifi- 


cant deviation from equilibrium 


The sulphur compounds the light oil were 
thiophenic character. The ratio 3-methylthio- 
phene the toluene fractions was examined. 


this work typical high sulphur heavy Canadian 
crude from the Lloydminster field was chosen 
feed material, and detailed study was made the 
gaseous and liquid products pyrolysis from inert 
fluidized solids system. The effect reactor temperature 
and residence time the composition these phases 
was particular interest. 

The significance this study lies 
production, from low cost raw material, gas rich 
unsaturated hydrocarbons and liquid product contain- 
ing large proportions aromatic compounds. 

Previous studies Peterson have shown that 
substantial upgrading similar heavy crudes may 
obtained thermal cracking under mild conditions (500- 
550°C.) fluid-solids reactor. Conventional cracking 
catalyst gives better results under such conditions but 
suffers from heavy coke deposition and metal deactiva 
tion. High temperature pyrolysis variety oil frac 
tions has also estigated Reid ‘al (2) using 
fluid bed, yields natural gas substitute were 
obtained this means. The large scale application the 
fluid bed technique for cracking heavy feeds has been 


confined the widely described fluid coking process 
(8. 4, 5, 6, 7, 8 


Gaseous product distribution resulting from the pyro 
lysis various hydrocarbon feeds has 


recently Linden and Peck who used empty 
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tubular reactor for their investigation. While such 
reactor permits establishment required residence times 
with considerably accuracy and between wide limits, 
does not allow accurate estimation true gas temper- 
ature, major variable which can measured directly 
fluid-solids system. 


Gas flow rates through fluid bed cannot varied 
greatly; however, possible induce reasonable 
range residence times varying bed depth and 
reactant recycle ratio. 


Experimental 
(a) Equipment: 

The continuous pilot plant used for this study has 
been described previously Figure simplified 
flow diagram the equipment. The preheated oil was 
injected directly into the bottom the reactor vessel; 
the usual practice feeding the reactor riser was not 
followed because rapid coke formation and consequent 
blocking the small riser. The properties 
minster crude oil are given Table The bed material 
used for all runs was finely divided silica sand; the screen 
analysis given Table 


PROPERTIES LLOYDMINSTER CRUDE 


Gravity 0.98 
Degrees 12.4 
Water Content w/w 0.2 
Conradson Carbon 8.2 


Distillation: 186-56 


20°, 14 
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TANK 


GAS HOLDER 


Figure 1—Pilot plant layout. 


Reactor temperature was automatically maintained 
supply hot sand from the burner vessel, which was 
operated temperature 50-100° higher than that 
the reactor. similar way differential pressure 
recorder/controller actuating valve the reactor 
standpipe, controlled the level the reactor vessel. Part 
the metered reactor offgas was recycled provide 
solids transfer the reactor and ensure smooth fluid- 
ization. 

addition product collection from the reactor, 
analysis samples were withdrawn intervals from the 
metered burner offgas and from the circulating burnt 
and coked sand stream. 


Operating conditions are summarized Table Al- 
though most runs lasted ten hours significance was 
attached analysis products collected during the first 
four hours. The range useable reactor bed depths was 
limited the fluidizing quality the material; very low 
beds created temperature control problems and high beds 
caused excessive slugging with attendant coke lump for- 
mation the reactor. The residence times given Table 
are average values calculated from the bed depth, bed 
void volume and from the mean reactor offgas flow rate 


TABLE 
SCREEN ANALYSIS BED MATERIAL 


Mesh Percent 


Weight 


Through 100 


Through 100 200 


Through 200 325... 


Through 


the temperature the bed. Use the latter parameter 
implies instantaneous generation the final gas volume 
the feed inlet, valid assumption the primary cracking 
reactions take place during small fraction the total 
residence time. Linden and Peck have made similar 
assumption. The recycling reactor offgas increased the 
single-pass residence time factor 1/1-R where was 
the fraction recycled. The walls the reactor above the 
fluid bed level were kept temperature 100°C. below 
that within the bed, minimize additional cracking out- 
side the reactor zone. 


(b) Results: 

Table summarizes the overall yields gas, oil and 
coke from each run. will noted that the weight yield 
process gas reaches maximum temperature 
about 650°C., while volume yields rise throughout the 
temperature range with formation larger proportions 
hydrogen and coke. 

The scatter data does not permit any conclusions 

reached about the effect residence time, nor 
does the process gas analysis, shown Figure reflect 
any trend this regard. The variations reactant partial 
pressure, caused changing bleed gas (nitrogen) pro- 
portion were not large enough effect greatly the 
general proportion the components the process gas. 

The total sulphur the light oil product decreased 
slightly with increase reaction temperature, and the 
same time there was corresponding increase the yield 
heavy sulphur residuum. The amount sulphur 
detected the process gas (as H,S) and the coke (as 
SO, the burner offgas) was but small proportion 
the total yield. 

Some typical A.S.T.M. distillation curves for the oil 
product are shown, with that the crude, Figure 
While the products from runs higher temperatures 
contain larger proportions light material, reference 
Table will show that such gains are partly offset 
lower overall yields. 
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Figure 3—Representative light oil distillations (ASTM 86- 
56) reactor temperature °C. 


(c) 
(i) Gaseous product distribution 


grams the corresponding values for thermodynamic 


equilibrium are plotted. The latter were calculated from 
librium” certain discreet reaction systems with resi- published free energies formation apparent 
times empty tubular reactor; marked simi- Figure that within the accuracy the data, 
wide variety feed materials. was interest see The reaction times involved have 
whether similar results were obtained the pyrolysis been sufficiently long being 
the components this system, and this respect, the 
Distribution reactants the following system was large positive deviations from found Linden 


considered: 


Figures and show reactant ratios for Systems (1) 
and (II) function temperature; the same dia- 
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TABLE 


Reactor 


Burner 


Bed 
Depth 
(inch 


Resid. 
Time 
(sec. 


Gas 


600 49 
650 
650 
650 36 
650 47 
675 
700 
700 
700 
725 
750 
800 


Recycle 


Bleed 
Rate 


Bed Air 
Depth Rate 
(inch) 


Bleed 
Gas 
.c.f.m. 


Offgas 


Offgas 


0.47 
0.37 
07 0.41 
0.37 
00 34 


404 
350 0 
374 
391 0.7 670 
675 
362 715 
750 
720 
710 
730 
740 
810 
750 
770 
820 
850 


640 
630 
675 


NN NNN NNN NNN 
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and Peck with residence times below about 1.5 seconds 
have not been duplicated. Some unpublished data obtained 
South American (Bachaquero) crude fluidized coke 
bed, fit remarkably well Figure and thus tend 
confirm that the behavior this system substantially 
independe ent feed material and the nature the 
inert bed used. 

The increased scatter data for System (Figure 
largely due the inaccuracy chromatographic 
estimation the small quantities propane the gas 
mixture. most cases there significant positive devia- 
tion from not apparent that increased resi- 
dence time would reduce this deviation, but since some 
the values depart almost order magnitude 
from median condition quantitative information 
obtained from Figure 

System has been con- 
perhaps very reaction. The results obtained 
this study confirm System was all 
very far from equilibrium. 
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(11) Liquid product composition 

The light oil products from several representative runs 
were closely fractionated end point 200°C. 

ft. mm. Podbielniak column atmospheric pres- 
sure. The overhead fractions were analysed for hydro- 
carbon types the method D1319); 
the data are summarized Figures and 

The steady increase aromatic content the gasoline 
fractions with increasing reactor temperature quite 
marked, reappearance small amount olefins 
800°C. has not been satisfactorily explained. Two other 
runs 650°C. were also fractionated and analysed see 
whether any marked differences resulted from variations 
reactor residence time; change type distribution 
was apparent and therefore the data are not shown. 


The light oil contained moderately large amounts (of 
the order thiophenic sulphur, the discreet frac- 
tions aromatic type were accompanied the thio- 
phene isomer corresponding structure. The sulphur 
content each fraction runs was 
determined; the results are shown Figure 10. Since 
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Figure sulphur light oil fractions. 


certain gaseous systems (viz. and can serve “tags” 
indicate roly sis severity, was thought that some 
analogous system the liquid phase might serve similar 


Accordingly, the proportions the two isomers, 


2-and were determined analysis 


the toluene cut from each fractionation using liquid 
partition chromatography and the ratio these iso 
mers, denoted f,, compared with the ratio expected 
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Figure 7—Distillation light oil fraction run No. 
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Figure 9—Distillation light oil fraction run No. 16. 
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REACTOR TEMPERATURE 


ratio with equilibrium 
values. 


shows that all samples tested there appears 
general trend for the ratio approach 
further supported the increase the ratios 
with longer residence times. However, apparent change 
reaction pattern has occurred the region near 
and thus more detailed interpretations are hardly war 
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Nomenclature 


absolute pressure, atmospheres. 
reaction temperature centigrade. 
partial pressure product 
partial pressure product 
mole fraction product, 
equilibrium. 
equilibrium. 
equilibrium. 


w/w percent weight. 


standard cubic feet per minute. 
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Empirical Correlation Azeotrope Data 


Extension the Meissner-Greenfield 


The Meissner-Greenfield correlations azeo- 
trope composition and boiling point have been 
checked with new data with satisfactory results 
except the cases hydrocarbons—acids and 
halogen hydrocarbons—acids. 


modification these correlations for azeotrope 
composition has been extended other combinations 
organic liquids where the deviations from Raoult’s 
law can predicted from the hydrogen bonding 
characteristics the components. 


simple empirical linear equation for azeotrope 
point has been applied most the cases 
included the above correlation. 


attempts have been made correlate 
with the thermodynamic properties solutions. 
general the more successful correlations are very empi- 
rical nature. 


Hildebrand has given general rules, that mini- 
mum boiling azeotrope possible when the two com- 
ponents have nearly equal vapor pressures and also when 
the two components differ considerably polarity even 
the boiling points are appreciably different and that 
maximum boiling azeotrope possible when there are 
negative deviations from Raoult’s law and small differ- 
ence boiling points. Kireev has pointed out that the 
smaller difference between the vapor pressures the 
two components given temperature the boiling 
range the more likely the existence constant 
boiling mixture. 

Lecat Rossini and Skolnik have pre- 
sented correlating equations and graphical methods for 
correlating azeotrope data where one component forms 
ith homologous series other components. 
These methods are quite accurate but are sometimes 
restricted. 

Meissner and Greenfield have proposed empirical 
correlations for azeotrope composition and temperature 
which are not dependent homologous series. Their 
correlation applicable systems here one component 
hydrocarbon halogenated hydrocarbon and the 
other alcohol, and aldehyde, ketone, carboxylic 
acid phenol. Their method has been extended this 
paper. 


Manuscript received October 30, 1958. 
‘Associate Professor, Department of Chemical Engineering, University of 
Toronto, Toronto, Ont. 
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Contribution from the Departinent of Chemical Engineering, University 
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Classification solutions according hydrogen 
bonding 


Hydrogen bonding refers the coordination 
drogen bonded carbon other elements 
with atoms such oxygen, nitrogen fluorine. 
Deviations solutions from ideality and hence azeotrope 
formation has been attributed large extent hydro- 
gen bonding. Ewell have classified liquids and 
solutions formed mixing these liquids the manner 
summarized Table This table indicates the classi- 
fications solutions covered Meissner and Greenfield 
their correlation and those which have been included 
this study. 

this paper the main source azeotropic data was 
the compilation Horsley 


Test the Meissner-Greenfield correlation 


Meissner and Greenfield proposed their correlation 

tions for azeotrope data compiled Horsley for the 
period 1947-52. Table lists the compounds included 
Figures 

these figures the azeotrope composition, 
mole fraction hydrocarbon and 


(1) 
where Boiling point hydrocarbon, 
Figure line was recommended Meissner for 
drocarbons and halogenated hydrocarbons with alco- 
hols only while line applied hydrocarbons with 
aldehy des, phenols, cresols, ketones 
and hologenated hydrocarbons with phenols and cresols. 
Figure are shown Meissner’s two curves for cor- 
relating experimental data hydrocarbons 
genated hydrocarbons with acids. 


For predicting azeotropic boiling point Meissner 
recommended the correlation the group 


3 
B 

with the difference boiling point between the high 
boiler and the low boiler the azeotrope com- 
bination. Here the boiling point the azeotrope. 
The value depended the system; typical values 
were tabulated Meissner. Figure shows the applica- 
tion their correlating line with new data. 


169, 
955). 
ons”, 
). 
edge- 
type 
irried 
iples, 


TABLE 
CLASSIFICATION ORGANIC ACCORDING HYDROGEN BONDING POTENTIALITY AND DEVIATION FROM LAw 


Class 


Class Characteristics Examples 


P-N 


Three dimensional network 
strong bonds 


Alcohols, acids, phenols, primary and second- 
ary amines, nitro-compounds with active 
atom, etc. ! 


| 


Both donor atoms and active 
atoms the same molecule 


Ethers, ketones, aldehydes, esters, nitro- 
compounds without active H-atoms, tertiary 
amines etc. 


Only donor atoms 


Only active atoms 


Halogenated hydrocarbons like 
etc. 


bond forming capabilities 


Hydrocarbons, sulphides, mercaptans, halo- 
genated hydrocarbons not class 


LEGEND positive deviations. 

negative deviations. limited solubility. 

(1) solutions covered modified Meissner correlation extended 
this study. 

(2) solutions covered Meissner and 


Halogenated 
Hydrocarbon 


Mole fraction 


0.2 


Halogenated 
hydrocarbons 


Recommended line for 

hydrocarbons and 

genated hydrocarbons 
with alcohols 


Recommended line for 
hydrocarbons with 
aldehydes, phenols, 
and cresols 


0.2 0.4 0.8 
Mole fraction hydrocarbon 


Mole fraction acid 
Figure Meissner’s correlation for 


and halogenated hydrocarbons with and for hydro- 

with aldehydes, phenols, and cresols with more 

boiling point hydrocarbon, and T,, boiling point 
other solvent, 


Figure 2—Test Meissner’s correlation for hydrocarbons 

and halogenated hydrocarbons with organic acids with 

more recent data. Here (T, where 

point other solvent, °K. 
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hydrocarbon 
Hydrocarbons 
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Table 


HALOG: HYDROCARBONS 
Dichloromethane 

Carbon tetrachloride 

Bromo methane 

bromo methane 

Bromoform 

Trichloroethylene 
Tetrachloroethylene 
Tetrachloroethane 
UNSATURATED Pentachloroethane 

2,3 Dimethyl 1,3 butadiene Cis bromo 1,2 dichloroethylene 
Heptyne Trans bromo 2,2 dichloroethylene 
3-Chloropropene 

1,2 Dichloro-1-propene 

Bromopropane 


DROCARBONS 


PARAFFINS 

Hexane 

2,3 Dimethyl butane 

Heptane 

Octane 

2,5 Dimethyl-hexane 
1,2,4 Trimethyl pentane 
Tridecane 


ALICYCLICS 

Cyclopentane 

Methyl cyclopentane 

Cyclohexane 2-lodopropane 

cyclohexane 1-Chlorobutane 

1,3 Dimethyl 
1-Bromobutane 

UNSATURATED 

Methyl cyclopentene 

1,3 Cyclohexadiene 

1,4 Cyclohexadiene 


AROMATICS ALCOHOLS 
Benzene Methanol 
Toluene Ethanol 


Butyl benzene alcohol 
Mesitylene Allyl alcohol 
Styrene Butyl alcohol 
o-oxylene butyl alcohol 
m-xylene Isobutyl alcohol 
Sec. alcohol 
alcohol 

amyl alcohol 
Iso alcohol 
2-Pentanol 


a-Terpinene 
Terpinolene 


PHENOS-CRESOLS ACIDS 

Phenol Formic 

o-cresol Acetic 

3,4 xylenol Propionic 
Butyric 

KETONES Valeric 
Isovaleric 


2-Butanone 
3-Pentanone 
4-Methyl-2-pentanone 


Extension Meissner-Greenfield correlation other 
H-bonding classes organic liquids 

For the combinations Table other than those 
treated Meissner and Greenfield, the Meissner-Green- 
field correlation was first attempted. some cases slight 
modifications the term were used. These modified 
forms are denoted y,, and The proposed 
equations were obtained the method least squares. 

Figure attempt has been made show the 
correlations obtained for systems which one com- 
ponent from class halogenated hydrocarbon. 

curve data for halogenated hydrocarbons with 
hydrocarbons were best correlated with Meissner’s 


with the proposed correlating equation 


curve data for halogenated hydrocarbons with 
halogenated hydrocarbons were correlated using 


obtain the correlating equation 


Curve shows limited amount data for halo- 
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genated hydrocarbons with water correlated using 
Meissner’s 


and curve shows small amount data with glycol 
correlated 


Curve shows large number points correlated 
using 


the correlating equation 


These points are for halogenated hydrocarbons with 
esters, ketones and on. 

Figure attempt has been made correlate data 
for class involving ketones and esters 
monocarboxylic acids. general 


Figure 3—Test Meissner’s azeotrope boiling point cor- 


3408 
and T,, and are the boiling points 
the low and high boiling pure components and the 
azeotrope, °K. 


O O 
C) 
is 


lowing simple expression can used. 
The proposed correlating equation 


Figure shows data for class V-V-hydrocarbons and 
drocarbons. Using the defined Equation (11), 
the data were correlated using 


Finally Figure shows data for class When the 
second component was water 2-aminoethanol the 
was used. With glycol and glycerol the 
following was used 


The correlating equation 


Table lists the compounds included Figures 3-7. 


Curve Classes 


Mole fraction halogenated hydrocarbon 


Figure 4—Azeotrope composition correlations for systems 

which one component halogenated hydrocarbon. 

where and T,, are the normal boiling points the pure 

components, °K. 


Correlation for predicting azeotropic boiling point 

The following simple linear equation was found 
boiling point data for both 
systems and most other drogen bonding classes 
which the azeotropic composition correlation was 
extended: 


and are constants which have been determined for 
each combination empirically and their values are given 
Table The best straight lines obtained 
squares from literature data are shown Figure 


Discussion 


The Meissner-Greenfield equations were checked 
using new data compiled Horsley and were found 
give agreement except for solutions hydrocarbons 
and hologenated hydrocarbons and acids. 

These correlations were extended, sometimes 
slightly modified form. The new correlations gave best 
results cases such halogenated hydrocarbons-donor 
atom liquids ketones-esters (although only few azeo- 
tropes hydrocarbons-hydrocarbons, and hydro- 
carbons-associated liquids the first class. The average 
deviation comparing actual azeotrope composition 
atmospheric pressure with the correlating equations 
(mole). New correlations have been indicated for 
other classes but the results were not good the 
above four classes. 

Data for systems not included Table showed 
tendency correlate the Meissner-Greenfield manner. 

azeotrope boiling point correlation has been pro- 
posed which simpler and applicable wider range 
compounds than Meissner’s. However not quite 
accurate Meissner’s method. 


0.2 0.4 0.6 0.8 1.0 
Mole fraction ketone 


Figure composition correlation for systems 
involving ketones and esters acids. Here 


the boiling points the pure components, °K. 
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ast 


While the graphs and equations presented permit the 
satisfactory correlation large amount azeotrope 
dat: they not predict the formation azeotrope; 
non azeotropes are reported Horsley for which 
azeotropes would predicted Meissner and the 


extension this paper. 


CLASS 
Water 
2-Aminoethanol 
Glycol 


Glycerol 


CLASS III 

KETONES 

Acetone 

2-Butanone 
2-Pentanone 
3-Pentanone 
3-Methyl-2-butanone 
2,4 Pentanedione 
2-Hexanone 
3-Hexanone 

4-Methyl-2 pentanone 
Cyclohexanone 
4-Heptanone 

5-Methyl-2 hexanone 
Acetophenone 

2,6 Dimethyl-4 heptanone 
Methyl-acetophenone 
Propiophenone 

Carvone 

Pulegone 


ETHERS 

Ethyl ether 

Propyl ether 
ether 
Ethyl propyl ether 
Butyl ether 

Benzyl ether 


ESTERS 

Ethyl formate 
formate 
formate 
formate 
formate 
Methyl acetate 
Ethyl acetate 
2-Methoxy ethyl acetate 
2-Ethoxy ethyl acetate 
Propyl acetate 
Isopropyl acetate 
Butyl acetate 

acetate 
Bornyl acetate 
propionate 
Ethyl propionate 
Butyl propionate 
Isobutyl propionate 
Methyl butyrate 
Methyl isobutyrate 
isobutyrate 
Methyl isovalerate 
Propyl isovalerate 
Methyl caproate 
Ethyl caproate 
lactate 
Propyl benzoate 
Ethyl a-toluate 


Table 


ALDEHYDES 
2-Furaldehyde 


ACETALS 

Methylal 

Acetal 

Acetaldehyde acetal 


OXIDES 
Propylene oxide 
oxide 


CLASS 
HALOGENATED HyDROCARBONS 
Iodomethane 
Dichloromethane 
Bromo-dichloro-methane 
1,2 Dibromomethane 
Chloroform 

Bromoform 

Carbon tetrachloride 

Bromoethane 

1,1 Dichloroethane 

1,2 Dichloroethane 
Trichloroethylene 

1,1,2,2 Tetrachloroethane 
Pentachloroethane 
1-Chloropropane 
2-Chloropropane 
2-Bromopropane 
1,2-Dichloropropane 
1,2,3-Trichloropropane 
1,3 Dibromopropane 

Chlorobutane 

Bromobutane 
1-Bromo-2 propane 
lodobutane 

butane 


CLASS 

PARAFFINS 

Hexane 

butane 

2-3 Dimethyl butane 
n-Heptane 
pentane 
2-3 Dimethyl pentane 
2-4 pentane 
2,2,3 Trimethyl butane 
2,2,4 Trimethyl pentane 
2,3,4 pentane 
2,5 Dimethyl hexane 

heptane 
Nonane 

Decane 

2,7 


PARAFFINS 
Cyclopentane 

Cyclohexane 

Methyl cyclopentane 

Ethyl cyclopentane 

Cis 1,1,3 trimethyl cyclopentane 
Trans 

Cis 1,3 Dimethyl cyclohexane 
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OLEFINS 
2,3,4 trimethyl 2-pentene 


CYCLO OLEFINS 


Cyclohexane 


AROMATICS 

Benzene 

Toluene 

Cumene 

Mesitylene 

Propyl benzene 
Cymene 
Pseudocumene 
Butyl benzene 
Styrene 

Ethyl benzene 
p-xylene 
Naphthalene 
naphthalene 
methane 
Acenaphthene 

1,2 ethane 


100 


-20 


TERPENES 
a-Pinene 
a-Terpinene 
nopinene 
Camphene 
Terpinolene 
Dipentene 
2-Limonene 


HALOGENATED AROMATICS 
Bromobenzene 
p-Dibromobenzene 
o-Dichlorobenzene 
m-Bromotoluene 

1-Bromo naphthalene 
1-Chloro-naphthalene 


Carbon disulphide 
Ethane thiol 

Allyl sulphide 

sulphide 

Butyl sulphide 

sulphide 


Mole fraction hydrocarbon 


Figure 6—Azeotrope composition correlation for solutions 
(T, where and are the boiling 


points the pure components, °K. 
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Figure 7—Azeotrope composition correlation for solutions 

hydrocarbons with Type solvents. Here (T, 

except for glycol and glycerol systems where 

components, °K. 


Nomenclature 


Intercept azeotropic boiling point correlation lines. 

Difference boiling point two components. 

equation constant (Graph 3). 

Boiling point, degrees Kelvin. 

Meissner’s term for correlating azeotropic boiling point. 

Mole fraction. 

Meissner’s term for correlating azeotropic Composition. 

Difference boiling points high low boiler from 
azeotropic boiling point. 


Subscripts 


High boiler 


hydrocarbon 
Halogenated hydrocarbon 
donor liquid 
Halogenated hydrocarbon 
hydrocarbon 
phenols 
ketone 
Hydrocarbon— acids 
Halogenated hydrocarbon-alcohols 
Hydrocarbon alcohols 
Halogenated hydrocarbon- water 
Hydrocarbon formic acid 


Figure 8—Correlation azeotrope boiling points. Here 
the difference between the boiling points the pure 
components and defined 

for minimum boiling azeotrope 
and T,, for maximum boiling azeotrope. 


TABLE 


Systems Classes 


Hydrocarbons formic acid. 
Hydrocarbons phenols V-II 
Hydrocarbons ketones 
Hydrocarbons V-V 
Hydrocarbons strongly associated 
Halogenated hydrocarbons donor liquids. 
Halogenated hydrocarbons 

hydrocarbons, 


Low boiler 
Azeotrope 
Modifications Meissner’s term. 
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The use solar energy, most applications, 
dependent primarily economic considerations; 
solar energy can become important energy source 
when the costs its use become competitive with the 
costs energy from other sources. The basic and 
usually most costly item equipment for solar energy 
utilization the collector, solar heat exchanger, 
and the amortization the first cost this and 
associated equipment the major cost solar 
energy. this paper, the probable economic feasi- 
bility several solar applications assessed 
using the costs energy from present sources 
for estimating upper limits the first cost 
solar heat exchangers and their associated equip- 
ment. For example, for space heating temperate 
climates, heat exchanger costs would have 
the range $1-6 for solar heating 
this appears within attainable range. Eco- 
nomic considerations have helped shape the develop- 
ment the Wisconsin solar energy program, and 
studies solar heat exchangers and their application 
heating, cooling and power generation are briefly 
described. 


OLAR radiation reaching the earth constitutes enor- 
mous energy resource. continuing energy 
supply which will available after the supplies readily 
available fossil fuels have become expensive 
cient for our needs. The solar energy supply several 
orders magnitude larger than the world’s total energy 
demands. The challenge solar energy research 
make the use this resource technically and economically 
feasible. 

The amount solar radiation incident the earth 
great, but the energy flux received the earth’s surface 
both low and variable with time and location. The 
intensity radiation varies from zero maximum 
about hp./sq. yd. and sunny, temperate climates, 
the annual average radiation horizontal surface 
the range 1500 2000 B.t.u./sq. ft./d. Table shows 
average daily June, December, and annual radiation for 


1Manuscript received October 23, 1958. 

Director, Solar Energy Laboratory, University of Wisconsin, Madison, 
Wis. 

8University of Wisconsin and Chemical Engineering Consultant, Denver, 
Col. 

Based on a paper presented at the Joint A..Ch.E.-C.1.C. Chemical Engin- 
eering Conference, Montreal, Que., April 20-23, 1958. 


Solar Energy 


Economics and Engineering Research 
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TABLE 


AVERAGE RADIATION FOR 
SEVERAL 


Annual 
San Juan, (18° Lat.) 1540 2000 1940 
Paso, Tex. (32° Lat.) 1210 2730 2040 
Fresno, Calif. (37° Lat.) 610 2620 1660 
Denver, Colo. Lat.) 670 2180 1550 
Madison, Wis. (43° Lat.) 450 1990 1220 


several locations. The collection and utilization this 
energy for most purposes are not now economically 
attractive, any solar-operated device system hich 
requires much energy must depend large areas 
energy collecting surfaces which are present expensive 
construct. The major problems solar energy utiliza- 
tion can considered technical and economic 
nature, although sociological other considerations may 
equally more important some applications. Dis- 
cussions these problems are given the proceedings 
several recent symposia 4), 

The sun radiates effective temperature about 
6000°C.; the radiation reaching the earth’s surface (after 
scattering and absorption some wave lengths 
atmospheric components) divided about the 
visible range, 51% the near infra-red, and the 
The radiation reaching the earth 


_consists diffuse (scattered) and direct (beam) com- 


ponents; cloudy days all the incident radiation may 
diffuse, and very clear days the total 
incident radiation may direct. The orientation 
surface has substantial the 
amount radiation intercepts, with surface contin- 
uously oriented toward the sun receiving the maximum 
radiation. temperate climates, the optimum orientation 
fixed receiver for winter-time use south facing, 
with tilt equal the latitude plus degrees, 
for summer use the optimum tilt less 
degrees. 


hols 
er 

ure 
ux. 

36, 


The most important part the solar energy utiliza- 
tion system the solar heat exchanger, “collector”. 
are two types, shown schematically Figure 

The flat-plate type uses solar energy absorbing sur- 
equal area the total area the collector, 
now apparently the less expensive build and can absorb 
both direct and diffuse radiation. The focusing type uses 
optical system, usually reflector, concentrate the 
direct component the radiation onto heat exchanger 
smaller than the projected area the collector, thus 
achieving higher energy flux the exchanger and reduced 
area for thermal energy loss. The flat-plate heat exchanger 
should have its most practical applications low temper- 
ature operations whereas the focusing collector should 
best applied where higher temperatures are required. 


principle, any process which utilizes heat could have 
this energy supplied from solar energy. However, the 
intermittent and diffuse nature solar energy make its 
utilization for some purposes more promising than for 
others. Included among those processes which appear 
have potential economic interest are comfort heating and 
air conditioning, power for special applications 
pumping irrigation water, and evaporation distillation 
saline waters recover potable water. these appli- 
cations, methods are available for handling the fluctuations 
the energy supply. The intermittency solar energy 
supply can ‘be “smoothed out” adding thermal inertia 
storage the comfort heating process, and good corre- 
lation between energy supply and cooling demand mini- 
mizes the storage requirements for space cooling. 
applications such pumping irrigation water solar 
distillation, the equipment can sized produce the 
desired total output during periods operation, and 
inertia then has the form product storage. Other 
applications potential interest are solar cooking and 
food refrigeration. 

Two other types processes interest are photo- 
electric and photochemical processes. The photoelectric 
processes permit direct conversion radiation elec- 
trical energy semiconductor cells, described 
Pearson Photochemical processes which utilize energy 
the solar spectrum endothermic reactions that can 
subsequently reversed recover the energy are 
interest the collection and storage solar energy. 


FOCUSING FLAT- PLATE 
COLLECTOR COLLECTOR 


TRANSPARENT. 


COVERS 


SURFACE 


REFLECTOR 


CONDUIT 
HEAT TRANSPORT FLUID 


solar collectors. The focusing collector operates only 

direct, beam, radiation. The flat plate collector operates 

direct and diffuse radiation, and does not have 
moveable. 


Examples this type reaction, although still the 
research stage, are the photochemical decomposition 
water the presence cerium iron salts obtain 
hydrogen and oxygen, studied Heidt and photo- 
chemical decomposition nitrosyl chloride, studied 
Niewirth and Daniels 


ECONOMICS SOLAR ENERGY UTILIZATION 


Solar energy most applications not now com- 
petitive with energy from other sources, and until is, 
widespread use can not expected. The costs collect- 
ing and using solar energy include the major items 
amortization first cost, direct operating costs, and 
realistic maintenance costs. There have been little quan- 
titative data available these costs for most the 
promising applications. Economic studies can serve 
indicate the more promising applications, where research 
and development studies may fruitful. 


The first cost system for utilizing solar energy, 
for any the purposes indicated, includes (a) the col- 
lector system with its pumps, ducts, piping, etc., (b) the 
storage system, required, (c) the “energy-utilizing” 
unit and (d) the necessary auxiliaries such control 
systems and instruments. These first costs should 
amortized over periods which realistically represent the 
expected life the equipment. Direct operating costs 
should include power costs for such purposes pumping 
heat transport fluids and labor for operation. Maintenance 
expense includes costs cleaning and repair collector 
systems, glass reflector cleaning and replacement, ser- 
vicing plumbing duct-work, and maintenance 
storage units, engines, coolers, and auxiliaries. 


Where continuous output continuous availability 
output process required, most temperate climatic 
zones appears economically impractical provide 
sufficient collector area and storage capacity operate 
processes exclusively solar energy. Auxiliary energy 
sources must provided carry the load during pro- 
longed non-sunny periods. ith most processes, the 
major economic consideration the minimization the 
total cost the product the process. Thus, fuel 
the major cost factor conventional system and the 
costs the product -operated system are 
lower than the corresponding fuel costs the conven- 
tional system, total costs can reduced use the 
solar-operated system. For example, house heating 
system areas where central heating required, 
conventional furnace needed; the total cost heat 
delivered from solar heating system lower than the 
fuel saved it, the overall cost heating reduced 
use the solar heating system. 


Costs energy for several applications 


The availability and costs energy 
sources over the world are widely variable, 
possibilities for economic utilization solar energy are 
correspondingly variable. The competitive 
solar radiation energy resource for many applica- 
tions may expected improve with time, advances 
the new fields solar technology make its use less 
expensive and the costs energy from other sources 
rise. The economics solar energy utilization are also 
dependent the development local resources. The 
competitive position solar energy thus varies with 
location, nature the application, and time, and general- 
izations the relative costs energy from solar radi- 
ation and other sources are difficult make. 
discussion, the probable economic feasibility several 
solar applications assessed using the range costs 
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Location Type and Capacity Cost 


TABLE 


Costs ENERGY FOR SEVERAL APPLICATIONS SELECTED LOCATIONS 


Notes and References 


HEATING 


Central U.S. 105B.t.u./hr. oil burner 


Western U.S. LPG system 


Natural gas 


South 2-5 tons compression 


South 2-5 tons absorption 


Cairo, Egypt tons compression 


Dakar, West Africa tons compression 


POWER 


Industrial large capacity 


Small central power plant 


Rural U.S. Auxiliary generator 


Cairo, Egypt 
Dakar, West Africa 


Rural India Animal power 


energy from present sources basis for estimating 
upper limits the first costs the solar operated 
systems. 

Table shows the range costs energy for three 
applications which are potential interest 
processes. 

Solar and conventional processes are compared 
estimating the maximum cost chargeable solar- 
operated process order for compete with con- 
ventional energy sources. The comparison simplified 
two assumptions. First, even though solar intensity 
widely variable over the earth’s surface, the climates 
where solar energy most likely widely used 
not have greatly differing annual solar energy reception. 
Second, although various designs and applications involve 
solar collection efficiencies which vary considerably, re- 
search materials and thermal properties collectors 
should lead improvements design and selection 
systems which permit average efficiencies conversion 
solar radiation sensible latent heat orking 
fluid reach about 60%. This figure used the 
following analysis, except where another choice dictated 
the special nature the application. Attainment 
this efficiency contingent further engineering ad- 

vances collector design and operation. 

solar energy applications discussed here, only the 
fixed costs, including amortization the solar equipment, 
are considered. Since there information available 
the literature maintenance operating costs for any 
solar operated systems, they have not been included. 
Ultimately these costs must minimized satisfactory 
design. 
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Fuel oil 15¢/gal., fuel cost only. 
Propane 17¢/gal., fuel cost only. 


Fuel cost only. 


Electric power 20% load factor, 
equipment cost included 10%/yr. 


Fuel cost only, $0.75 
Same, with equipment cost included No. 


Electric power 7.5¢/k.w.h. fixed cost. load factor. 


Electric Power 17¢/k.w.h. 
factor. 


vr. fixed cost, load 


Power generated large central thermal power stations 
hydroelectric stations. 


Power from 10,000 k.w. central power station. 


Power from small gasoline generator plants 
capacity. 


(7) 


Residential lighting rate 
Residential rate 


Power for lifting irrigation water estimated terms 
electrical equivalent. 


Some the assumptions made estimating the allow- 
able costs solar units are part arbitrary, and the 
estimates are considerably affected thereby; this partic- 
ularly true the assumptions useful life the equip- 
ment. Where there reasonable indication the useful 
service life solar unit, this has been taken the basis 
for amortization cost, the amortization rate has been 
taken equal pical rate for the energy source 
for investment the solar operated systems. 


contrast many items process equipment, the 
cost per unit capacity for system flat plate solar 
produced units, would nearly independent the 

capacity the system. Thus, the estimates followi ing, 
the costs are considered the basis square foot 
collector without regard the total capacity the 
system. The assumption that cost independent 
capacity for focusing collector less valid because 
variability requirements for moving and positioning the 
reflector, wind loading, and other factors. 


Space heating 


Space heating application solar energy which 
interest because (a) the temperatures required are 
low enough (100° 150°F.) permit the use flat-plate 
collectors, (b) the solar energy incident dwelling 
during the winter most temperate climates more than 
the energy needed for comfort heating, (c) there 
requirement for conversion thermal other forms 
energy, and (d) comfort heating major user 
conventional fuels. 
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Space heating costs the United States are large 
extent dependent the proximity the user fuel 
supplies. Table the approximate costs one million 
B.t.u. net heat delivery residence from inexpensive 
natural gas (based rates Denver, Colo.), propane 
“bottled gas” (typical western U.S. price) and fuel oil 
national price) are listed. estimating space 
heating costs, combustion efficiency 80% has been 
assumed, and only fuel costs have been considered the 
fixed costs the conventional equipment are 
not reduced the solar heat supply. the United States, 
typical costs fuel for residential space heating generally 
range from 75¢ $2.00 per million B.t.u. delivered. 
other areas, costs may several times great. 


With heat collection one square foot flat-plate 
solar heat exchanger rate 150,000-300,000 B.t.u. per 
annual heating season (variable with the length the 
heating season), the value the output per square foot 
exchanger would range from 11-60¢, based value 
delivered heat $0.75-$2.00 per million B.t.u. fixed 
cost expenditure $1.10-6.00/sq. ft. could 
justified for the collector and associated equipment. 
The flat-plate heat exchangers suitable for house heating 
applications can fabricated from materials such 
glass, sheet metal, and insulation; costs these materials 
most proposed designs total less than $1.00/sq. ft. 
collector. However, the labor costs for fabrication the 
solar collectors for the few solar houses built date have 
been many times the material costs. 

For practical solar heating system, heat storage must 
also provided. This may the form liquid 
(hot water), solid (crushed rock), change-of-state 
material. With heat storage capacity equivalent 
day’s solar collector delivery capacity, the cost heat 
storage equipment should much less than that the 
solar collector. Auxiliary equipment such fans, ducts, 


piping, and control instruments, where required addi- 
tion the conventional system will add further invest- 
ment the system. However, appears that collectors 
and associated equipment can built, quantity, for 
costs within the range $1.00-$6.00/sq. ft., and develop- 
ment solar space heating systems appears attractive. 


Space cooling 

Research solar space cooling suggested the 
generally good correlation between energy supply and 
cooling demand, areas high solar incidence usually 
have uncomfortable summer temperatures. Technical 
feasibility assured the fact 
absorption refrigeration units require temperatures well 
within the range working temperatures solar heat 
exchangers. Comfort cooling solar energy 
associated with solar heating system, and used for 
both purposes, the allowable cost the solar heat ex- 
changer would greater than for either the separate 
purposes because its useful energy collection over 
greater part the year. 

evaluating the potential economics solar cooling, 
costs energy for conventional refrigeration systems 
system has been evaluated terms total fixed and 
operating cost with use factor 20%. For absorption 
systems, which heat could supplied the refrigerant 
generator either solar conventional means, 
interest make two comparisons: the fuel cost only can 
compared with the cost heat from solar collector, 
the total costs both fuel and equipment 
compared. costs for air conditioning vary 
greatly with location shown the examples selected, 


which were chosen primarily for their 
high solar availability, and also their wide energy cost 
range. 

Table shows typical Unites States costs the range 
refrigeration, applicable small com- 
pression absorption units two five tons refrigera- 
tion capacity. heat-operated system the cost fuel 
would usually lie the range equivalent 
$36 $72 worth fuel per year 1200 operating hours. 
solar heat substituted for fuel absorption 
system about 150 sq. ft. solar collector would 
necessary provide sufficent heat for replacing each ton 
cooling capacity. Considering use solar energy only 
for cooling, 10% annual fixed cost, $2.50 $5.00 could 
paid per square foot solar collector and the neces- 
sary storage facilities. These costs are approachable 
range. Substitution solar energy for the usual heat 
source conventional cooler will require some modi- 
fication cooler design, particularly flat-plate collectors 
are used, because the variable nature the solar energy 
supply. 

Replacement entire compression cooling system 
solar operated absorption type cooler requires con- 
sideration air conditioner costs well solar heat 
exchanger costs. total conventional costs 5-8¢/t.hr. 
refrigeration, there would investment approxi- 
mately $4.00-$7.00 allowable per square foot solar col- 
lector and associated air conditioning equipment. the 
costs air conditioning unit operable with solar heat 
are reducible those conventional unit, $100-$200/t. 
refrigerated capacity, the allowable cost the solar 
collector would reduced the range $3.00-$6.00/sq. 
ft. 

Two other factors improve the prospects 
cooling. many installations, particularly temperate 
climates, space heating will also accomplished use 
the solar collector; this increases substantial 
amount the allowable investment the solar collector, 
storage unit, and other common facilities. Secondly, in. 
other parts the world where energy more expensive 
than the United States, the economic situation more 
attractive for the solar operated cooling system. Dakar, 
for example, the cost electricity for air conditioner 
operation ten times great the United States. 
This permits substantial increase allowable investment 
solar cooling facilities. many such areas space cooling 
now used very limited extent. the purchasing 
power the peoples these countries rises, however, 
and the modern conveniences enjoyed industrialized 
countries come into increasing demand, solar 
ditioning should achieve substantial application. 


Solar power 

Generation power from solar energy, the form 
mechanical work, has been accomplished limited 
scale the production steam other vapor 
boiler heated solar collector and the expansion 
this fluid reciprocating engine. Electricity has been 
generated coupling generator this type engine, 
converting sunlight directly electrical energy 
silicon semi-conductor cell “solar battery”. These 
applications have been either experimental specialized, 
and not based the economics substantial energy 
supply. 

Costs power from several selected types power 
installations shown Table range 
Most residential consumers electricity the United 
States and Canada pay for power from central 
plants. Power from some small central installations 
more expensive, ranging and that from 
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diesel-powered generator sets used non- 
farms for pumping irrigation water, etc., costs 
range 4-8¢/k.w.h. some foreign locations 
much higher costs prevail, even central stations. 
these costs, was necessary some cases 
assuine load factors and fuel costs; the best available data 
and reasonable assumptions load factors and useful life 
were employed. most the small instal- 
fuel the major item cost. 


solar power plant consisting focusing 
collector flat- solar producing steam 
for use expansion engine, the maximum plausible 
efficiency from steam (100 p.s.i.) electricity small 
units appears about 15%. With focusing collector 
having 60% solar collection efficiency, about sq. ft. 
solar collection surface would required for the 
daily generation one k.w.h. Over the total cost range 
electric power shown Table from 1-30¢, and 
total annual fixed cost 15% the original invest- 
ment, the allowable expenditure for solar collector, power 
plant equipment, and all auxiliary facilities would range 
from $24 $720 for the sq. collector and 
associated equipment necessary develop one 
Per square foot solar receiving surface, these figures 
become $1.60 and $48, respectively; for average 
delivery one kilowatt for hours daily, $240 $7200 
are the maximum allowable costs for the 150 ft. 
collector and associated power plant. 

Small steam engines, manufactured quantity, 
should obtainable for less than $50/k.w. capacity, and 
electric generators should the $50 range. $100 
for these items, the allowable investment solar heat 
exchange units would $140 $7100. Thus, areas 
high power costs (e.g., rural India and parts North 
Africa), the allowable investment solar boiler and 
associated equipment high enough make solar power 
potential economic interest. With these high allowable 
investments per square foot solar reflector flat-plate 
exchanger, the prospects appear however, 
there are small prime movers presently available that 
will operate efficiency near 15% 100 steam, 
nor are there focusing collectors reasonable cost which 
will produce this steam 60% efficiency. further 
deterrent these areas may lack capital make 
the initial investments required for the solar power 
systems. 

Since solar power not continuous 
storage appears prohibitively costly, auxiliary 
energy supply would have provided “where con- 
tinuous power essential. The load factor the auxiliary 
facilities would reduced the solar power plant, 
and effect, only portion the fuel cost reduction 
the competing (standby) electric power station should 
credited the solar power plant. some applications, 
such pumping irrigation water, auxiliary energy ould 
not required and the economic comparison involves 
the total alternative original investments and operating 
costs. 

For large central power plants, has made 
economic studies the use flat-plate collectors and 
low pressure turbines and has concluded that the costs are 
potentially the range two three times the current 
costs power from conventional fuels the United 
States. 

Other processes 

There are several additional uses solar energy which 
merit consideration. Among these are applications non- 
industrialized areas where. sociological considerations may 
important the economic factors; solar cookers and 
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food coolers are examples. solar cookers can made 
cheaply enough and the requirements for their opera- 
tion can made consistent with the eating and living 
habits the users, they could used millions 
people who are now using wood other primitive fuels. 
northern Mexico, for example, fuel for cooking 
major item expenditure for many low-income families, 
and reduction the use wood could improve their 
economic situation well decrease the rate de- 
forestation. Small, inexpensive refrigerators, operated 
solar energy other heat sources, could provide facilities 
for food preservation not now existent many these 
areas. 


Salt water demineralization solar distillation could 
provide potable water arid regions where solar radia- 
tion high. distilled water cost $2.00/1000 gals. 
(based costs several large fuel-operated evaporation 
plants) and typical daily productivity one-tenth 
gallon potable water per square foot solar distillation 
about cents the maximum allowance for 
investment each square foot solar distiller and 
accompanying auxiliary equipment annual fixed and 
operating costs are 10%. With the lower cost capital 
used for public works, annual fixed costs could some- 
what lower and the total allowable investment might 
increased $1.00-$1.25/sq. ft. Large scale construction 
methods and the use cheap materials such bituminous 
coatings, glass, and some plastic films may permit con- 
struction this cost range. The problem difficult, 
however, and exceptional economies design, construc- 
tion, and operation must developed solar distillation 
compete with fuel-operated demineralizers. 

comparison with the usual costs natural fresh 
water, the projected prices for solar-distilled water are 
very high. However, the need for demineralized saline 
water expanding, and favorable locations, solar 
distillation saline water has the potential becoming 
economically attractive where potable water 
available for human consumption. 


Summary 

The foregoing appraisal several solar energy appli- 
cations from the standpoint their competitive position 
with respect conventional sources energy shows 
the possibilities and limitations these applications. 
each case chosen, solar energy appears have promise 
some areas and some scales use. Other applications 
may similarly considered. 


THE WISCONSIN PROGRAM 

The Engineering Experiment Station the University 
Wisconsin, cooperation with the chemical, mechan- 
ical, civil engineering, Chemistry and meteorology de- 
partments, studying some the problems the 
engineering applications solar energy. This program 
includes studies flat-plate solar heat exchangers and 
focusing collectors, intermittent 
absorption cooling systems, focusing collectors for small 
power units and food cookers, photochemical processes, 
and solar radiation measurements and correlation. 

The Wisconsin program was initiated 1953, when 
conference persons who had been active research 
this field was held Madison. This conference was 
supported the National Science Foundation and the 
University. meeting, and result the 
interest Professor Farrington Daniels the Chemistry 
Department, program research the field was 
initiated. With support from the John Simon Guggen- 
heim Foundation, the National Science Foundation, the 
Wisconsin Alumni Research Foundation, and the Uni- 
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versity, and more recently with grants from the 
Foundation, several principal investigations, and 
number related researches engineering, 
and meteorology have been undertaken 


Collector research 


Basic all the thermal applications discussed the 
solar heat exchange unit, either the flat-plate 
focusing type. The economic studies outlined indicate 
that the major item cost space heating and cooling 
and small power systems the collector. Performance 
this unit described heat balance for the solar 


energy absorbing surface, which can expressed as: 


HR = Qu + Cth Che 


where 


solar radiation rate horizontal surface (total 
radiation for flat plate heat exchanger, direct radiation 
only for focusing 
the ratio radiation tilted surface that 
horizontal surface. 
projected area the collector. 
the solar energy absorbing surface area the solar 
heat exchanger for flat-plate exchanger). 
effective specular reflectivity the reflector 
focusing collector. (This factor does not enter for 
flat-plate exchanger). 
the transmissivity for solar radiation any system 
transparent covers over the absorbing surface. 
absorptivity the blackened surface for solar radiation. 
the useful energy transferred the heat transport fluid 
the collector, per unit area absorber, per unit 
tuume. 
the sum all losses from the absorber, including losses 
convection and radiation through the transparent 
covers, and downward losses conduction, and shad- 
ing, per unit area absorber, per unit time. 

heat capacity effect the heat exchanger, per unit 
area absorber, per unit time, not operating 
constant temperature. 


The left side the equation represents the absorbed 
radiation, while the terms the right show the distribu- 
tion the absorbed radiation. This heat balance has been 
extensively studied for flat-plate heat exchanger 
Hottel and Woertz For either type collec 
performance improved minimizing the optical losses 
the departure the product from unity) and 
the thermal losses, the most 
usually the losses from the absorber convection and 
radiation. 

The objective collector studies minimize the 
cost collected energy, achieving higher collector 
and/or use more economical materials 
and designs. The Wisconsin flat-plate heat exchanger 
studies are presently oriented toward the reduction 
the optical losses, which limit the 
solar energy, and the reduction upw ard radiation losses 
from the absorber the use semi-transparent coatings 
make surfaces having selective radiation characteristics. 
Progress has been made the latter area, and surfaces 
oxides metals, having absorptivities for solar 
radiation above 0.9 and emissivities for long wave radia- 
tion less than 0.4, are now being used 
experimental solar heat exchangers. The recent 
ability plastic materials with desirable characteristics 
plate exchanger covers for fabrication 
collector parts has permitted study their performance 
and utility. The concept unitized solar heat exchanger 
design, suitable for factory assembly and easy installation, 
being explored means reducing costs. 


Focusing collector studies have included consideration 
the use plastics and plastic forming techniques for 
transparent coatings applicable focusing collector 
receivers also, they are designed operate tem- 
perature range above 250°F. where radiation losses may 
become large fraction the thermal losses. Several 
geometric arrangements reflector and receiver are 
being studied, are the precision fabrication and 
orientation the reflectors (which limit the ratio 
areas reflector receiver). The overall objective 
minimize the cost collected energy the temper- 
ature levels required for promising applications. 

Figure shows experimental flat-plate exchanger 
under construction, and Figure shows experimental 
focusing collector with flat-plate receiver. The perfor- 
mance both types has been discussed 
publication 

Solar cooling 

The Wisconsin program solar cooling development 
includes consideration continuous absorption cycles 
food coolers. Intermittent units, designed for 
food refrigerators non-industrialized countries, have 
been described Williams, The continuous 
systems under study are absorption coolers with heat 
supplied the generator the cooling units 
plate solar heat exchangers. schematic diagram shown 
Figure storage unit shown between the collector 
and the generator; alternative the location storage 
facilities between the aporator and the cooled space. 


assembly, prior blackening the copper absorber sur- 

face. Provision made for mounting three glass covers 
over the absorber. 
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for with polystyrene-aluminum reflector and flat plate heat 
exchanger. 

Ambient temperature, horizontal radiation, and radia- 

tion incident south- -facing tilted collector during 

typical hot summer day Madison show Figure 
the variability radiation fixed flat-plate solar heat 

exchanger and approximation the changing cooling 

loads. study the relationships betw een cooling load 

and energy supply collectors various orientations 

indicates that hot days Madison are usually accom- 

panied high radiation. 

The challenges this program lie the development 
techniques meet the limitations and variability 

the solar energy supply suitable design the com- 


solar exchanger- cooler system. The immediate 
objectives are investigate some the many possible 
design combinations solar heat exchangers, short-term 
storage systems, and cooling units, analy tical studies 
considerable importance optimize collector design for 
combined summer cooling and winter heating use. 


Solar power 


The development small h.p.) solar power 
units using steam engines other prime movers being 


CONTINUOUS SOLAR ABSORPTION COOLER 


GENERATOR 


SOLAR 
COLLECTOR, 


COOL AIR 
ROOM 
STORAGE 


Figure diagram continuous solar absorp- 
tion cooler. 
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studied. The problems are divided into separate investi- 
gations the solar exchangers and prime movers, with 
subsequent combination the two elements obtain 
optimum systems. The solar heat exchanger studies out- 
lined above contribute the solar power investigations, 
and parallel with them have been studies several con- 
ventional and specialized types prime movers which 
might combined with ‘solar collector. Most these 
efforts engines have been concerned with low pressure 
steam, p.s.i., and higher pressure ranges are now 
being investigated. The experimental results with low 
pressure engines indicate that with presently available 
equipment this capacity, engine efficiencies are less than 
4%; 15% will very difficult achieve unless sub- 
stantially higher steam pressures 
better solar heat exchangers) are employed. For solar 
power units underdeveloped areas, ease maintenance 
and operation unskilled personnel also being con- 
sidered the development and design equipment. 

Another approach the problem solar power 
generation lies the use solar “batteries” light- 
sensitive silicon cells. The present cost these units 
prohibitive, but the temperature the cells can 
controlled, there the possibility increasing the energy 
flux them use focusing reflector. With the 
increased electrical output per cell, the number cells 
required and the total investment may reduced. Studies 
these factors are progress. 


Other applications 


applications solar energy under- 
developed areas the economic factors may more 
important considerations. This point 
may illustrated reference recently developed 
reflective type comprising conveniently 
supported plastic parabolic reflector large enough 
deliver heat grill rate 600 800 watts. Some 
these units have been distributed low-income fam- 
ilies northern Mexico with the objective assessing 
their effectiveness replacing part the wood 
other cooking fuels used that area. 

connection with problems obtaining solar radia- 
tion data the desired form, i.e., the total radiant 
energy over hourly periods, integrator-recorder has 
been dev eloped. This instrument shows hourly values 
integrated radiation rather than the record radiation 
time usually obtained with recording potentiometers. 
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Figure 5—Ambient temperature, radiation horizontal 
collector, and radiation collector tilted 35° the 
south, for hot summer day Madison, Wis. 
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Conclusion 


The challenges solar energy research are great, and 
the potential fruits engineering research and develop- 
ment studies this field are substantial, they can 
achieve the major objective making solar processes 
economically competitive with processes using other 
energy resources. Solar energy tremendous, virtually 
untapped resource which should available long 
men are earth use it. The Wisconsin program seeks 
contribute the solution the problems which must 
met before this energy resource can significant 
factor our economy. 
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The Effects Shape and Reynolds 
Number Drag the Motion 


Freely Oriented Body Infinite Fluid’ 


The drag freely oriented bodies motion 
relative infinite fluid related the Reynolds 
number, the formula 


where Re, formed from the diameter, 
D,, sphere with the same surface area the 
ber. The inertial drag coefficient, C,, related 
shape, orientation, and the Reynolds number. 
Reynolds numbers above 2000 the relation ex- 
pressed the formula, 


where the surface sphericity, the ratio the 
surface area sphere with the same volume 
the surface area the body. Correlations are also 
developed for predicting Reynolds numbers 


below 2000. 


motion relative infinitely extending fluid has 
been solved analytically only few special cases. 
Dimensional analysis was therefore used Rayleigh 
obtain the well known drag formula, 


The unique relation between the drag coefficient, f(Re), 
and the Reynolds number, Re, predicted this formula 
for geometrically similar bodies similarly oriented with 
respect the motion has been investigated for wide 
variety bodies. However, few attempts have been 
made correlate the effects shape and orientation 
they affect chemical engineering operations. Since this 
problem importance operations such sedi- 
mentation and elutriation and techniques such sedi- 
mentation analysis, was decided subject 
general study. 
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Theory 


The Rayleigh drag formula, Equation (1), states the 
general predictions dimensional analysis concerning 
the solution the drag problem. does not, however, 
incorporate explicitly any expectation the form 
the relation between the drag coefficient and the Reynolds 
number. Since such expectation exists, (as shown the 
following evident that this general statement the 
drag law capable some refinement. 

The Navier-Stokes equations—the general differential 
equations motion viscous fluid—contain viscous 
and inertia terms describe the features real fluid 
motion, and expected that the solutions 
these equations are also composed additive viscous and 
inertia terms. This shown, for example, Oseen’s 
approximation their solution for sphere, which 
inertia term added Stokes’ familiar solution for 
the viscous drag. The additivity viscous and inertia 
drags also suggested principle the Rayleigh drag 
formula itself, for region where the inertia forces 
are unimportant the drag must proportional the 
viscous forces, while region where the inertia forces 
dominate the reverse approximately true. may there- 
fore supposed that the drag body motion 
through infinite fluid composed viscous drag 
and inertial drag: 

The proportionality between the viscous drag the 
viscous forces formally stated the Stokes’ -type 
expression 


~ 


which defined the viscous drag coefficent, 
(the coefficient the viscous drag force), and 
arbitrary. Similarly, proportionality the inertial 
drag and the inertia forces gives Newton’s formula 

where the drag coefficent. inserting 
these values the component drags Equation (2) and 
multiplying through there obtained the 
quadratic drag formula 


which defined the drag number. may noted 
that Rayleigh’s total drag coefficient, (f(Re) Equa- 
tion(1)), and Reynolds number are based the 
same characteristic diameter then and 
f(Re).Re?. also noted that the viscous 
drag coefficient definition independent the Rey- 
nolds number, since any drag excess that predicted 
Equation (3) directly ascribable the action 
the inertia forces. The viscous and inertial 
cients are coefficients terms drag equation, but 
should not confused with Rayleigh’s total drag co- 
efficient defined Equation (1). The only case 
which equality exists high Reynolds numbers, when 
the inertial drag coefficient becomes formally equal 
Ray leigh’ drag coefficient the dimen- 
sion expressed consistent terms. 


may seen that the above quadratic formulation 
the drag law has the following advantages over Equa- 
tion (1): 

The coefficient constant for given shape and 
orientation, and the coefficient should only weakly 
dependent the Reynolds number. contrast, the 
total drag coefficient f(Re) strong function the 
Reynolds number. 

The equation gives quantitative picture the 
viscous and inertial aspects drag. This likely 
first importance the transition zone, and may facilitate 
rational inv estigation the effects shape orien- 
tation this region. 

The drag and the fluid velocity not appear 
together any term, whereas the conventional drag 
coefficient involves Moreover, the dimensions 
the body appear only the Reynolds number, 
terms involving the drag force. 


Different definitions the charcteristic dimension 
appropriate the nature the drag, can used the 
Reynolds numbers, (right side Equation (5)). Simi- 
larly, the coefficients and can related different 
parameters describing shape and orientation. This ad- 
vantage will most apparent the transition zone. 


The remaining problem, before applying the quadratic 
formula the analysis experimental data decide 
suitable parameters for describing the size and shape 
solid bodies motion through fluid. evident 
that, order have practical merit, these parameters 
should both simple and universal. This criterion 
met reasonably well parameters constructed from the 
volume the surface area and the area projected 
the direction the motion. There are thus obtained 
(i) the surface sphericity 


and the form sphericity 


where For prolate bodies there may also 
added the ratio, (the length diameter ratio). 
Previous investigators have attempted relate the 
total drag coefficient Equation one another 
these There not appear any other 
criteria which are worthy practical consideration when 
dealing with freely oriented bodies. 


may observed, conclusion this section, that 
general the inertial drag, (Equation (4)), can re- 
solved into component resulting from tangential forces, 


Fir = V6Ciz pulo*As........ 
and component resulting from normal forces, 


where For bluff shaped bodies tangential 
(friction) drag generally small compared with the 
normal (form) drag, while for streamlined bodies the 
opposite true. Most chemical engineering problems treat 
with bluff shaped bodies, and hence this distinction not 
present interest. 


The viscous drag coefficient 


Stokes, neglecting the inertia terms the equations 
motion viscous fluid, derived theoretically the 
law, verified numerous experimental measurements 
(5.2, 8), for the viscous drag sphere: 


The value, (in Equation (3)), the viscous drag coeffi- 
cient for sphere is, thus, 24. 


the general case, which analytical solution 
the equations motion not practicable, the viscous 
drag body arbitrary shape can expected 
proportional the intensity shear the surface and 
the surface area. If, then, the diameter the body 
taken the diameter, sphere with the same 
surface area, the variation the viscous drag coefficient 
with shape and orientation should small and should 
directly related the attitudes the surface the 
motion. will therefore assumed that 


where the function equal unity for spheres. 

Figure shows for oblate and prolate bodies, 
evaluated from the data Heiss and Coull and for 
isometric bodies, evaluated from the data Pettyjohn 
and Christiansen The oblate and prolate bodies were 
oriented, respectively, with their flat sides long axes 
either parallel perpendicular the motion. The iso- 
metric bodies showed effect orientation 
terminal velocities still fluid and give values 
unity can whence 


practically interesting consequence this result 
that the volume and density small particle and its 
terminal velocity gravitational field are known, 
Equation (12) can used estimate its surface area. 


The inertial drag coefficient 


The inertia drag sphere very low Reynolds 
numbers has been approximated theoretically 
who obtained for the inertial drag 
value 


SPHEROIDS 

CYLINDERS 

>» RECTANGULAR PARALLELOPIPEDS 
ISOMETRIC POLYHEDRONS 


Figure 1—The shape function for oblate, prolate and 
isometric bodies. 
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Figure 2—Inertial drag coefficients freely oriented isometric bodies and disks function the Reynolds number. 


high Reynolds numbers, the inertia forces pre- 
dominate over the viscous forces everywhere except 
the shallow boundary layer and the wake. Here 
Prandtl’s boundary lay theory has been used estimate 
the tangential drag variety bodies, chiefly 
streamline shape (5), However, this theory has little utility 
treating the bluff, freely oriented bodies present 
interest: the flow around bluff body normally separates, 
and the greater part the drag due the normal 
stresses produced this separation. 


The inertial drag sharp-edged, bluff body 
separating flow can expected roughly propor- 
tional the area projected the direction the motion. 
However, chemical engineering operations 
tional complexity introduced the absence orien- 
rotational and other motions the general translation 
the body relative the fluid. The projected area 
which the inertial drag freely rotating body should 
roughly proportional not clearly defined, but 
appears reasonable choose the projected area, Ag, 
sphere equal surface area; this area roughly inter- 
mediate between that sphere equivalent volume 
and that the sphere swept out free rotation the 
body. inserting this value the projected area the 
expression for the drag, Equation and noting 
that the viscous drag best expressed putting 
Equation (3), the exact formulation (for the present 
purposes) the quadratic drag law, Equation (5), be- 
comes 


Figure shows inertial drag coefficients for isometric 
bodies, calculated Equation (13) from the data 
Pettyjohn and Christiansen and for disks, calculated 
from the data Squires and Squires function 
the Reynolds number. Reynolds numbers greater than 
five the inertial drag coefficients these bodies fall 
between the solution Blasius for flat plate oriented 
parallel the motion and the experimental value for 
disk oriented perpendicular the motion. Rey nolds 
numbers smaller than five the coefficient for spheres 
approaches Oseen’s solution, but the coefficients non- 
spherical bodies appear approach considerably smaller 
limiting values. The experimental evidence indicates that 
freely oriented bodies show three distinct types be- 
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havior, each which characteristic well defined 
range Reynolds numbers: 


10. Reynolds numbers smaller than 
the flow wholly laminar, and has been shown 
that bodies possessing three perpendicular axes sym- 
metry are equally stable all orientations The trend 
the behavior the inertial drag coefficient appears 
directly opposite that found all subsequent 
regimes; Figure indicates that the coefficient decreases 
value with increasing bluffness and angularity. 

(2) 100 400. Entry into this region 
marked the manifestation specific orientational 
stability the position maximum drag. Thus, isometric 
bodies and disks orient themselves with flat side forw ard, 
perpendicular the motion has been shown 
that Reynolds number about five the streamlines 
behind disk broaden out, and stationary region its 
rear separated from the main flow vortex sheet. 
this region retarded fluid there permanent 
circulation corresponding vortex ring. Hence, the 
appearance specific orientational stability 
bodies coincident with the development vorticity 
the wake, and directly associated with the stability 
the vortex ring. The dimensions the ring increase with 
increasing Reynolds number until critical Reynolds 
number reached which the stability the ring and 
orientational stability the body are lost together. For 
disks, critical Rey nolds number 195 has been found 
The orientational stabilities observed Pettyjohn 
and Christiansen indicate that the critical value may 
high 400 and low 100, but generallv 
the neighborhood 200. 


Figure indicates that the relation between the 
inertial drag coefficient and the Reynolds number this 
region accurately described the power formula 

Cy = Co. (14) 
The coefficients Figures and have common value 

Table gives the experimental values the exponent 
appears that the exponent lies between minimum 
value zero and maximum value 0.5, and pre- 
dicted with reasonable accuracy the formula 


1s not ov ec aa, a 8 a a a aa 
molds 
Oseen 
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Figure 3—Inertial drag coefficients freely oriented iso- 
metric bodies function the Reynolds number, 
showing the linearity the logarithmic relation the 
region where the orientation stable the position 
maximum drag. 


(3) 100 400. has been shown that when 
the Rey nolds number exceeds the critical value, approxi- 
mately 200, disturbance the vortex ring becomes 
visible, and successive portions its substance are dis- 
charged downstream regular intervals time. Freely 
oriented isometric bodies simultaneously develop rota- 
tional motion and notes that disks develop 
obbling teetering motion. The discharge vorticity 
increases rate and irregularity with increasing Rey nolds 
number, and eventually wake becomes completely 
turbulent mass eddies and vortices. Figure indicates 
that the discharge vorticity generally accompanied 
increase the inertial drag coefficient which con- 
tinues until the wake becomes completely turbulent 
Reynolds number near 2000. From this point onward the 
inertial drag coefficients sharp-edged bluff bodies are 
substantially independent the Reynolds number. The 
final portion the curve for spheres, Figure based 
the results Bacon and Reid and shows the 
decrease the coefficient the Reynolds number in- 
creases above about 105. This decrease typical 
round-edged bodies, and due the retardation 
separation induced the growth turbulence the 
boundary layer 


Experimental 
Pettyjohn and Christiansen defined the total drag 


and found that Reynolds numbers 


TABLE 
VALUES THE EXPONENT FOR FREELY ORIENTED BODIES 


Cube octahedrons 1.07 0.352 
Tetrahedrons 0.670 0.670 0.253 0.297 
Disks 0.000 0.000 0.000 0.000 


greater than 2000 the total drag coefficients 
freely oriented isometric bodies are related the surface 
sphericity the formula 


(18) 


when the viscous drag negligible, evident that 
purely inertial regime flow Also, appears 
from Equation (18) that 


(19) 


However, examination Equation (17) indicates that 
Equation (18) can only valid for high surface spheri- 
cities Equation (18) predicts maximum value 
5.31, whereas for flat, platy particles negligible 
thickness the factor D,? the denominator can 
become small that approaches infinity, indeed 
does Similarly, Equation (19) predicts value 
value. evident, therefore, that Pettyjohn and Christ- 
iansen’s formula cannot extrapolated beyond their data 
values smaller than those covered isometric 
bodies, 0.67 1.00). 

study the inertial drag coefficients non- 
isometric bodies was therefore undertaken test the 
generality Pettyjohn and Christiansen’s formula high 
surface sphericities and investigate the relation 


Shapes studied. Solid prisms and cylinders were cut 
from cork and cardboard. Table summarizes the pro- 
perties typical specimens. 


Method. The time fall stagnant air was determined 
from heights 840 and 1250 cm. 


TABLE 
DRAG COEFFICIENTS NON-ISOMETRIC BODIES FOR REYNOLDS 
Res 


Dimensions* 


Cylinders 


2.20 0.86 0.81 0.99 


0.231 0.35 0.94 1.82 


Rectangular Prisms 

0.180 1.00 0.80 0.89 1.34 


0.60 3.85 0.63 3.85 0.146 2.20 


Right Triangular Prisms 


2.50 7.6 0.063 0.172 0.95 1.82 


denotes diameter; width; height; length; and base. 
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Orientational behavior. The flat-shaped disks and 
when dropped, were oriented perpendicular 
the translation, and showed regular fluttering motion 
about the axis least inertia. The stability this orien- 
tation was destroyed atmospheric disturbances. 
The stable motion for all bodies consisted spiral 
translation accompanied rotation about the axis 
inertia, this axis being oriented perpendicular the 
translation. This motion could imparted the begin- 
ning fall manual rotation the body. The time 
fall flat shaped bodies with rotational motion was 
5-10% greater than for oriented fall the position 
maximum drag. 


Caleulation drag coefficients. The force balance 
for body falling from rest 


When the Reynolds number large, the drag can 
expressed 


The differential the motion is, therefore, 


Solution the differential equation, subject the initial 
the motion 


where 


whence, substituting Equation (23) for 
t 


tion the integral the right this equation then 
the approximation 


for values greater than five, (u, 0.985 
5). thus obtain for the drag coefficient 


The significance the drag coefficient calculated 
from the above equation follows. Until approaches 
value about 0.6, Equation (25) reduces 


(28) 


0.6, 0.3 u,. Hence, the terminal Reynolds 
number sufficiently large, can assumed that the 
Reynolds number will large when the point reached 
when the drag appreciably reduces the acceleration, and 
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the velocity this point will the same order the 
terminal velocity. The drag coefficient calculated from 
Equation (27) should, therefore, closely approximate the 
inertial drag coefficient for steady motion the time 
parameter and the terminal Reynolds number are both 
large. the present experiments was all cases larger 
than five, and the terminal Reynolds number was above 
2000. 


Table gives experimental values the inertial drag 
coefficient for number representative bodies. Figure 
shows the coefficient function the surface spheri- 
city for all the bodies investigated. seen that the data 
for high sphericities are fairly closely distributed about 
Pettyjohn and Christiansen’s data isometric bodies. 
significance can attached the unevenness the 
distribution, since not improbable that simply 
reflects the average nature nonisometricity among the 
shapes studied rather than general trend. The deviations 
from curve through the data isometric bodies are 
endurable from practical point view, and hence there 
justification for attempting recognize the effects 
nonisometricity quantitatively. Nor there usually 
any convenient basis for making such recognition. may 
concluded, therefore, that the present data are best 
represented curve which (1) begins the point 
0.42 0.44, 1.00 for the sphere, the point 
which all data must converge the surface sphericity 
approaches unity; (2) passes through the data iso- 
metric bodies; and (3) satisfactorily represents the data 
nonisometric bodies the range 0.67 hich lies 
below that 0.67) spanned isometric bodies. The 
simplest relationship which satisfies these conditions 


The present data indicate that the above relationship 
good for all freely oriented, bluff shaped bodies 
uniform density, with the exception pronouncedly 
prolate bodies; i.e., bodies with slimness (length dia- 
meter) ratios higher than about three. body with 
slimness ratio four was included the present studies, 
and gave value nearly 50% greater than predicted 
the equation. 

Flat plates and disks held perpendicular 
the direction fluid motion have inertial 
efficients which are practically identical with those 
freely oriented bodies given Equation (29); e.g., for 
square plate 2.32, while for disk 2.24, 


© PRESENT DATA 
@ ISOMETRIC BODIES 


Figure 4—Inertial drag coefficients isometric and non- 
isometric bodies function the surface sphericity 
Reynolds numbers greater than 2000. 
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INERTIAL DRAG CHARACTERISTICS FREELY ORIENTED BLUFF BODIES 


Flow regime Wake 
0.1 Viscous Irrotational 
0.1 5.5 Transition 


Transition 


5.5 200 Fixed vortices 


drag 


200 500 


vorticity 


500 3000 Transition Increasing disorder 


Fully turbulent 

tia, 
comes turbulent. 
Wakes rounded 
bodies narrow 


Inertial 


Inertial 


All orientations are stable when 
there are three more perpendi- 
cular axes symmetry 


Stable position maximum 5.5 


Transition Periodic discharge Unpredictable. Disks and plates 
tend wobble, while fuller bluff for 
bodies tend rotate 


Rotation about axis least iner- 
frequently 
spiral translation 


Free orientation 


compared with the value 2.25 predicted the 


equation. The slimness (length width) ratio rec- 
tangular plates has little effect until exceeds value 
ten. The inertial drag coefficient then climbs 
limiting value 3.9 slimness ratio infinity. These 
figures give some indication the behavior 


expected prolate bodies free orientation. 


may remarked that the total drag coefficients 
plates held perpendicular the fluid motion are con- 
ventionally based the projected area f(Re) 
and when the viscous forces are negligible, 
Thus f(Re) defined equals one half 
2f(Re) purely inertial regime and surface 
sphericity zero. 


Summary correlations 

The present paper relates the drag freely oriented 
body (or particle) motion through infinite fluid 
the relative fluid velocity and the properties the 
fluid and the body the relation 


| | | | 
4 + } 
0.2 


Inertial drag coefficient 


Inertial drag negligible 


For angular bodies assume 2.25. 
For spheres and rounded 
Figure 


For non-stagnant fluid assume 
2.53 when this gives the 
greater value 

0.8 


See Figure 


with 


Coefficients rounded bodies decrease 
the wake narrows 


where (It may noted that the com- 
monest value the drag force that obtaining when the 
body terminal velocity uniform gravitational 


field; Vg(ps-p))- 


The best experimental value the viscous drag co- 
efficient for randomly oriented bodies (excepting 
prolate bodies with slimness ratios larger than 
24. This value should accurate within +18% for 
individual bodies, and within few percent average 
for statistical populations. particular orientation 
individual body being considered accurate value 
junction with Figure 


Table gives summary formulae relating the 
inertial drag coefficients freely oriented bluff bodies 
shape and Reynolds number. The behavior the 
coefficients also shown graphically Figure The 
transition regime has been divided into four subregimes 
and the inertial regime into two subregimes. There 


| 
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Figure 5—General correlation the inertial drag coefficients freely oriented bodies with Reynolds number and shape. 
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disco..tinuity passing from transition subregime 
Goldstein notes that entry into the 
transition subregime III considerably advanced the 
turbulence the fluid, and for this reason 
stagnant fluids the correlation for the subregime III 
should extended subregime when this gives the 
higher value the coefficient. The correlations for the 
inertial drag coefficient not appear reliable for 
rolate bodies with slimness ratios greater than three. 
Other than this, they can applied all freely oriented 
bluff bodies. They should usually accurate within 
+20%, regions abrupt transition excepted (such that 
100 Reg 400). 


The correlations this paper apply motion 
infinite fluid. apply them motion cylindrical 
vessels may that the correlations Francis 
and Monroe for viscous and turbulent flow 
respectively are applicable bluff shaped bodies 
general: 

and 

where the ratio the diameter the 
body the diameter the column. 


The correlations this paper apply the motion 
isolated bodies. The quadratic drag formula, Equations 
(2-5), was first used the author Master’s thesis 
presented the University Saskatchewan the De- 
partment chemistry and chemical engineering, and 
entitled “Flow fluids through beds particles the 
fixed and fluidized states”, (May, 1955). The thesis (p. 
112) gives general correlation the relative terminal 
velocity settling ratio u/u, with the fraction void and 
the drag number for dispersions uniform particles 
motion relative liquid. The correlation facilitates 
extension the present results fractions void between 
0.45 0.55 and unity full spectrum Reynolds 


numbers. 


Conclusion 

The present work shows that the quadratic formu- 
lation the drag law facilitates precise, theoretically 
significant analysis experimental data 
oriented, bluff shaped bodies. Application the formula 
demonstrates the role the inertial forces low Rey- 
nolds numbers, and shows that data spheres indeed 
extrapolate exactly Oseen’s theoretical solution. 
interesting note that the coefficient the inertial drag, 
C,, decreases with increasing Reynolds number through 
the transition zone. This decrease connected with the 
growth the fixed vortices which attach themselves 
any bluff body this region. 

does not appear that the scheme correlation 
adopted the present study can simplified altered 
with any significant advantage; i.e. the coefficients the 
drag equation, presently defined, show minimum de- 
pendence shape and Reynolds number. The corre- 
lations extend over full spectrum shapes, orientations, 
and Reynolds numbers for freely oriented bluff shaped 


bodies, and thus provide practical information regions 
previously neglected, especially the transition zone. 
The weak dependence the coefficients and 
shape and Reynolds number, contrast previously 
proposed total drag coefficients, increases the reliability 
the present correlations and simplifies interpolation. 

The chief need for new supplementary data appears 
lie the following two areas: (1) data freely 
oriented prolate bodies turbulent flow; and (2) data 
non-isometric bodies the transition zone region 
with better expression for the exponent the 
Reynolds number. 


Nomenclature 


Dimensions absolute units: 


mass 
length 
time 
area projected body the direction its 
area projected sphere with the same surface 
viscous component the drag force.............ML/t? 
inertial component the drag 
local acceleration due gravity.............. 
velocity body relative the undisturbed 
fluid motion the body through infinite 
x/t average velocity the body over time 
terminal (steady-state) velocity approached 
the body after falling from L/t 
Dimensionless factors 
viscous drag coefficient 
inertial drag coefficient 
exponent dependent shape and orientation 
Rep 
Res 
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